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Abstract 


Purpose - Find optimal flight and design parameters for three objectives: minimum fuel con¬ 
sumption, Direct Operating Costs (DOC), and environmental impact of a passenger jet air¬ 
craft. 

Approach - Combining multiple models (this includes aerodynamics, specific fuel consump¬ 
tion, DOC, and equivalent CO 2 mass) into one generic model. In this combined model, each 
objective’s importance is determined by a weighting factor. Additionally, the possibility of 
further optimizing this model by altering an aircraft’s wing loading is analyzed. 

Findings - When optimizing for a compromise between economic and ecologic benefits, the 
general outcome is a reduction in cruise altitude and an unaltered cruise Mach number com¬ 
pared to common practice. Decreasing cruise speed would benefit the environmental impact 
but has a negative effect on seat-mile cost. An increase in wing loading could further optimize 
the general outcome. Albeit at the cost of a greater required landing distance, therefore limit¬ 
ing the operational opportunities of this aircraft. 

Research limitations - Most models use estimating equations based on first principles and 
statistical data. 

Practical implications - The optimal cruise altitude and speed for a specific objective can be 
approximated for any passenger jet aircraft. 

Social implications - By using a simple approach, the discussion of optimizing aircraft opens 
up to a level where everyone can participate. 

Value - To find a general answer on how to optimize aviation, operational and design-wise, 
by using a simple approach. 
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1 Introduction 

1.1 Motivation 

Optimizing aircraft and flight parameters for economic reasons goes back more than a centu¬ 
ry. People have been trying to minimize fuel consumption and the cost of flying since the in¬ 
vention of the very first aircraft. And nowadays the search still continues for more fuel- 
efficient engines, lighter structures or optimized routes for airlines. Economic reasons drive 
most, if not all, decisions in aircraft design and in operating procedures for airlines. 

In the most recent decades, following the growing worldwide attention to climate change, the 
amount of research on the environmental effects of aircraft has ramped up significantly. This 
is further augmented by the recent focus of regulatory organizations (e.g. ICAO) on decreas¬ 
ing the climate impact of the aviation industry. Nevertheless, attention to the environmental 
impact of aviation was and still is lagging behind on the automotive industry. 

The main issue with existing research is that optimization models are mostly kept behind 
closed doors by manufacturers and airlines to guard them from the competition. This makes it 
difficult for outsiders to gain any insight in realistic operating data of modern passenger jet 
aircraft. The aim here is to open up the discussion to more people by working with first prin¬ 
ciples, simple equations, and publicly available data. 

This thesis is written with a reader with a basic engineering background in mind. All aviation- 
related terms and abbreviations are defined in the list of definitions and list of abbreviations or 
elaborated in the body of the text. 


1.2 Title Terminology 

“Conditions for Passenger Aircraft Minimum Fuel Consumption, Direct Operating Costs and 
Environmental Impact ’’ 

Passenger aircraft: An aircraft which main goal is to transport passengers. In this thesis pas¬ 
senger aircraft will be narrowed down to passenger jet aircraft. 

Fuel Consumption: The mass of fuel burned per unit of time or distance. 

Direct Operating Costs: The costs directly associated with operating the aircraft. 


Environmental Impact: The direct or indirect impact on the environment of operating an air¬ 
craft. 
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1.3 Objectives 

The main objective of this thesis is to deliver an overview of the influence of different flight 
and design parameters on an aircraft’s fuel consumption, direct operating costs, and environ¬ 
mental impact. The influence of cruise speed Mach number, altitude, and wing loading will be 
analyzed. The aim is to make these models as generic as possible by minimizing the input of 
aircraft specific data. 

The first step is to review and analyze some fundamental methods for calculating fuel con¬ 
sumption, DOC, and environmental impact separately and the methods needed to calculate all 
relevant parameters. All models are based on first principles, simple equations, and statistical 
data. 

The models are then implemented in an Excel tool which calculates each of the previously 
mentioned goals and the most important parameters to reach them. The influence of each of 
these parameters is further analyzed and compared to real-world flight operation practices 
when possible. 

The final step is to combine all separate models into one model to find and optimum cruise 
Mach number and altitude for minimum fuel consumption, DOC and environmental impact 
combined. These results are then recalculated for the same aircraft with a higher wing loading 
in order to analyze the impact of this design parameter. The final result is a combined value 
which can be altered for greater economic or ecologic importance by means of weighting fac¬ 
tors. 


1.4 Previous Research 

Several cruise speed options to obtain minimum fuel consumption are discussed by Young 
2018. Specific air range is defined and calculated and by use of this parameter, three different 
cruise speed options are considered and compared with each other. Furthermore, an optimum 
cruise altitude for minimum fuel consumption is defined as cruise-climb and the commonly 
used stepped approximation of the cruise-climb maneuver is discussed. 

An important parameter for calculating fuel consumption is the thrust specific fuel consump¬ 
tion. The TSFC’s dependency on flight altitude, speed and thrust setting is discussed in Ben- 
sel 2018. A model for calculating the TSFC is made by Herrmann 2010 and is used in the 
fuel consumption model in later chapters. 

Direct operating costs have been researched by a multitude of organizations. For the model in 
this thesis, a method by Scholz 2015 will be used which is mostly based on the AEA1989 
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methods. This method is fully elaborated in the lecture notes by Scholz, therefore only results 
and some relevant equations will be mentioned in this thesis. 

The environmental impact of the aviation industry has gained more attention recently, provid¬ 
ing enhanced research. The first comprehensive report on this subject was written by the 
IPCC, called “Aviation and the Global Atmosphere” (Penner 1999). It handles the direct and 
indirect effects of the aviation industry on the climate. These effects are put into quantitative 
results by using radiative forcing, albeit with rather high levels of uncertainty. A differentia¬ 
tion is made between subsonic and supersonic aviation and their estimated future impacts. 
The report also lists multiple options to reduce the impact of aviation. 

One area where the IPCC report lacks in is accuracy of the impact of aircraft-induced cloudi¬ 
ness. This includes contrail fonnation and aircraft-induced cirrus clouds. In an update 
(Sausen 2005) to the IPCC report, the radiative forcing of aircraft-induced cloudiness was 
strongly reduced, while increasing the accuracy of the results. In a report by Fichter 2005, the 
influence of cruise altitude on contrail fonnation is discussed. 

The master thesis by Van Endert 2017 introduces the concept of ecolabels for aircraft. Alti¬ 
tude dependency of NO x and contrail fonnation are calculated and the concept of equivalent 
CO 2 mass is used as an alternative metric for radiative forcing. Most equations and data used 
were originally produced by Schwarz 2011, but have been simplified and compiled to gener¬ 
ate one value in the fonn of equivalent CO 2 mass. 


1.5 Structure of the Work 

Chapter 2: Overview of existing research and results as an introduction to the subject. 

Emphasis on fuel consumption and environmental impact. 

Chapter 3: Brief review of the models used in this thesis for aerodynamics, TSFC, 

DOC, the atmosphere, and equivalent CO 2 mass. 

Chapter 4: Explanation of the Excel tool used to calculate all results for Chapter 5 to 8. 

Chapter 5: Model for minimal fuel consumption and its most influential parameters. 

Chapter 6: Model for minimal direct operating costs and its most influential parame¬ 

ters. 


Chapter 7: Model for minimal environmental impact and its most influential parame¬ 

ters. Split into resource depletion and engine emissions. 
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Chapter 8: 


Chapter 9: 


Combined model of DOC and environmental impact. Includes a subchapter 
about the influence of wing loading. 


Conclusions and recommendations. 
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2 Literature Review 

2.1 Cruise Speed for Minimum Fuel Consumption 

The minimum fuel consumption case can be seen as a maximum range case. Minimum fuel 
consumption is considered as burning a minimum amount of fuel for a maximum distance, 
while maximum range can be described as covering a maximum distance on a minimum 
amount of fuel. These two definitions boil down to the same idea. 


For normal commercial aircraft operations, the cruise flight phase spans the majority of the 
total flight time. Taking this into consideration, most models only calculate optimum speed 
values for cruise flight. In addition, climb and descend speed can only be varied lightly due to 
different restricting parameters (e.g. ATC restrictions or runway length). 


Young 2018 gives an overview of different possible cruise speeds. These cruise speeds are 
based on the specific air range r a , which can be defined as: 


dx 

dx dt V V V 

—drrif —drrif W f c -T N c ■ D 

dt 


( 2 . 1 ) 


From this equation the theoretical speed for maximum r a , noted as V MSR , can be derived. The 
specific air range depends on the true airspeed V, the TSFC c and the drag D. The TSFC c is 
assumed constant and the drag D can be expressed in function of speed V: 


D = 


^DO ‘ P 'S 


V 2 + 


2 ■ (m ■ g) 2 
tt ■ A ■ e ■ p ■ S 


1 

Y l 


( 2 . 2 ) 


Or simplified: 


D = A D ■ V 2 + B 


D 


1 

Y 


(2.3) 


To find the highest SAR in function of speed (with a TSFC assumed constant), the ratio V/D 
should be maximized. The following derivation is made: 


d 

dV 



= 0 


(2.4) 


This results in: 
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V, 


MSR 


= V3- ^=V3-K 


md 


(2.5) 


With: F md the speed for minimum drag (or maximum aerodynamic efficiency). 


Continuing with Eq. (2.1), r a can be defined based on the TSFC. By integrating the resulting 
equation for r a the still air range R can be achieved. 


V ■ E 


c ■ m ■ g 


1 f ml V-E 
R — — ■ - dm 

g J m2 c-m 


( 2 . 6 ) 


(2.7) 


With: m 1 the start-of-cruise aircraft mass 

m 2 the end-of-cruise aircraft mass 


2.1.1 Analytical Integration 

An analytical integration of Eq. (2.7) can be done if the assumption is made that the TSFC 
does not change during cruise. The result is a reasonable approximation, but assuming a con¬ 
stant TSFC is the most important shortcoming in these simplified models for cruise speed and 
altitude. In the following integrations a mean constant value for the TSFC will be used, noted 
as c. 


The approach used by Young 2018 and Hull 2007 is to first establish boundary conditions for 
three different scenarios. Integrations can be found in the aforementioned literature, only the 
relevant results in terms of operation of the aircraft will be mentioned here: 


1) Cruise at constant H and constant C L 

Considering the equation for the lift coefficient in steady level flight: 

2 ■ m ■ g 
° L ~ p ■ V 2 ■ S 


( 2 . 8 ) 


It is clear that to keep the altitude (related to p ) and C L constant, V should decrease with 
decreasing m. Which means the engine thrust Fn must also decrease constantly, making 
this an unviable solution for the flight crew. 
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2) Cruise at constant V and constant C L 


This situation is known as cruise-climb. The resulting expression after integrating is 
known as the Breguet range equation: 


V-E 
R = -— 
c-g 



(2.9) 


For non-steady level flight Eq. (2.8) turns into: 


C L 


2 • m • g ■ cosy 
p-V 2 -S 


( 2 . 10 ) 


During cruise climb the flight angle y is very small, which leads to cos y = 1. As V and 
C L are to be held constant, the ratio m/p should also remain constant. This implies that 
p (and thus the altitude) should increase with decreasing m. Consequently, the altitude 
should be increased slowly during the whole cruise phase, a cruise-climb maneuver. 


3) Cruise at constant H and constant V 

In this case the integration can be solved if it is assumed that the drag can be modelled 
by the transonic drag polar (Eq. (3.1)). During cruise the lift coefficient in this scenario 
will reduce with time (due to its dependency on mass, which in turn decreases with the 
consumption of fuel). This means that the drag coefficient will also reduce. Considering 
the following equation for aircraft drag during cruise: 

T N = D=^-p-V 2 -S-C D (2.12) 

A reducing drag coefficient means reducing drag, which in tenn means the crew must 
constantly reduce the thrust to maintain a constant V. 


2.1.2 Maximum Range Cruise Speed 

Young 2018 describes the following three cruise speed options. The first is the maximum 
range cruise speed: MRC. This is the speed that results in the maximum SAR, which is (theo¬ 
retically) V MSR — V3 ■ V md . This speed decreases slightly while the aircraft becomes lighter 
due to fuel consumption during the flight, as illustrated in Figure 2.1: 
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For: 

• 35000 ft 

• ISA conditions 


Low gross 
weight 


0.65 0.70 0.75 0.S0 0.85 


Mach number, M 


Figure 2.1: MRC and LRC speed (Young 2018) 


The MRC speed depends on an aircraft’s gross weight, cruising altitude and the ambient air 
temperature and is the speed which corresponds with the lowest fuel cost. 


2.1.3 Economy Cruise Speed 

The cruise speed used in practice by most airlines is called economy cruise speed: ECON. It is 
a little faster than the MRC speed, therefore using a small amount more fuel but also arriving 
at the destination in a shorter time. Time is an important factor in calculating the trip cost as 
the time-dependent costs can be substantial. The ECON speed is usually calculated (by itera¬ 
tion) before each flight with the most recent cost index (Cl) from the airline (Roberson 
2007). This cost index is different for each operator and aircraft. 

The ECON speed is the speed for the lowest direct operating cost (and not just the lowest fuel 
cost). This time dependency of the direct operating cost can be seen in the seat-mile cost, 
modelled in Chapter 6. 


2.1.4 Long-Range Cruise Speed 

The long-range cruise speed LRC is a simplified approach to the ECON speed. For calculat¬ 
ing the LRC speed it is not necessary to have any knowledge of the airline’s cost structure. 
The LRC speed is determined by allowing a 1% reduction from the peak SAR and selecting 
the faster of the two resulting speeds (Airbus 2002). The advantage is that 1% of range is 
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traded for 3 to 5% higher cruise speed (Hurt 1965), which in turn leads to a shorter flight 
time. The difference between MRC and LRC speed is illustrated in Figure 2.2. 


2.1.5 Comparison 

The following figure illustrates how the different costs vary with increasing Mach number in 
cruise flight. As stated before, the MRC speed corresponds with the lowest possible fuel cost 
and the ECON speed with the lowest possible total cost. When time-dependent costs are con¬ 
sidered, ECON speed is always higher than MRC speed. The LRC speed is typically about 2- 
5% faster than the MRC speed according to Roberson 2007 and Seto 2009. 
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Figure 2.2: Cost overview for the different cruise speeds (Young 2018) 

Another thing to note is that the LRC speed is calculated for zero wind conditions, while the 
ECON speed usually takes wind speed into account when entered into the FMC. This means 
that the ECON speed will be lowered in a strong tailwind and increased when headwind is 
present, making optimal use of weather conditions (Roberson 2007). 


2.2 Cruise Altitude for Minimum Fuel Consumption 

The optimum cruise altitude is the altitude that results in the greatest SAR for a particular air¬ 
craft weight, ambient air temperature and Mach number. These data are determined by flight 
testing and published in the aircraft’s FCOM by the manufacturer. For many aircraft models, 
flight testing data shows that the optimum altitude tends to increase quasi linear with decreas¬ 
ing aircraft weight. To fly at the optimum altitude during the whole flight would mean in¬ 
creasing altitude constantly for the duration of the cruise (see Chapter 2.1.1, scenario 2). This 
maneuver is called cruise-climb and is forbidden by ATC due to flight level restrictions. 
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Flight levels (FL) are used to assure safe air traffic by separating aircraft into certain altitude 
zones. They are used starting from a certain transition altitude above the ground (depending 
on country). The altitude is calculated by measuring the ambient pressure and is standardized 
to the ISA pressure at sea level of 101 325 Pa. This standardization guarantees that every air¬ 
craft perceives flight levels at the same altitude, regardless of local variations in ambient pres¬ 
sure. Most modem aircraft are certified for reduced vertical separation minima (RVSM). This 
reduces the minimum flight level separation from 2 000 ft to 1 000 feet and allows aircraft to 
safely fly more optimum routes, reducing fuel consumption and increasing airspace capacity 
(FAA 2017). Even though the minimum separation has been reduced, aircraft are still not al¬ 
lowed to constantly change flight level during their flight. That is why in practice jet aircraft 
fly a stepped approximation of the cruise-climb maneuver, illustrated in Figure 2.3: 
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Figure 2.3: 


Step-climb approximation for cruise-climb (Young 2018) 


2.3 Direct Operating Costs 

The direct operating costs are ah costs made by an operator of an aircraft (e.g. airline) in order 
to operate the aircraft. The DOC can be split up into time-independent/accountancy costs like 
depreciation and interest, and into time-dependent costs. These time-dependent costs are 
called the cost index and include fuel, maintenance, crew, and fees and handling. Fuel price 
can be considered highly time-dependent (see Figure 2.4) and can alter a DOC model based 
on the date of calculation. 
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U.S. Gulf Coast Kerosene-Type Jet Fuel Spot Price FOB 
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Figure 2.4: Variation in jet fuel prices throughout the recent decades (EIA 2019) 


The method used in this thesis is described in the lecture notes by Scholz 2015. It is based on 
the methods defined in AEA 1989a and AEA 1989b. The AEA 1989a variation is for short- 
and medium-range aircraft and the b variation for long-range aircraft. The output can be rep¬ 
resented as the annual cost of one aircraft, denoted as: C a / ca or the seat-mile cost C s m . 


2.4 Environmental Impact 

The environmental impact of an aircraft in cruise flight can be divided into resource depletion 
and engine emissions with their impact on the global atmosphere. Noise pollution and local 
air quality will not be considered here since only cruise flight is covered in this thesis. 


2.4.1 Resource Depletion 

Resource depletion will be limited to the depletion caused by an aircraft’s fuel consumption. 
It could be seen in a broader context that includes the materials that were used to manufacture, 
transport and assemble different parts of an aircraft. However, the aim here is to treat the op¬ 
erational aspect of the aircraft only. Therefore, resource depletion can be set equal to the fuel 
consumption of the aircraft. Jet fuel (kerosene) is distilled from crude oil, which is considered 
the resource in this model. 


2.4.2 Aviation and the Global Atmosphere 

The environmental impact of aviation has been researched thoroughly by the IPCC in their 
1999 report titled “Aviation and the Global Atmosphere” (Penner 1999). The following 
summary is based on this report. 
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There are a number of different ways in which operating a jet aircraft impacts the environ¬ 
ment. Within the scope of this thesis only the impact of actually flying the aircraft will be 
considered, no indirect impact or emissions from ground transport at airports for example. 
First, there is the measurable emission from the combustion of jet fuel. This consists of car¬ 
bon dioxide (CO 2 ), water vapor (H 2 O), nitrogen oxides (NO x ) and other particulates like sul¬ 
fur oxides (SO x ) and soot. These gases and particles are emitted directly into the atmosphere 
around the tropopause level (see Chapter 3.4 for the tropopause). 

Subsequently, there is the climate impact of these emissions. This is a lot harder to measure, 
but by means of radiative forcing (RF) the climate impact of emission components can be 
compared mutually or to the climate impact of other industries. Radiative forcing (RF) is the 
difference between the energy which the Earth absorbs from the Sun and the energy which is 
radiated back into space. It is defined in units of watt per square meter. A positive radiative 
forcing will cause wanning of the Earth, a negative one will cause cooling. This balance be¬ 
tween absorbed and reflected energy determines the average global temperature and can be 
influenced by greenhouse gasses (among other factors). Radiative forcing makes it possible to 
measure climate change in a quantitative way. The drawback of using radiative forcing is that 
it cannot be applied to calculate the influence of a single flight, only the total influence of all 
aviation emissions during a period of time. In most cases this period is defined as 100 years 
(Jardine 2005). For this reason, the metric of equivalent CO 2 mass will be used (elaborated in 
Chapter 3.5). 

The following table gives an overview of each emission component’s direct or indirect impact 
on the environment: 
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Table 2.1: Climate impact of different emission components (Penner 1999) 


Emitted Species 

CO, 

Role and Major Effect at Earth's Surface 

Troposphere and Stratosphere 

Direct radiative forcing —* warming 

H,0 

Troposphere 

Direct radiative forcing —» warming 

Increased contrail formation —> radiative forcing —» wanning 
Stratosphere 

Direct radiative forcing —> warming 

Enhanced PSC formation —> O, depletion —» enhanced UV-B 

Modifies O, chemistry —> O, depletion —> enhanced UV-B 

NO x 

Troposphere 

O, formation in upper troposphere 
—»radiative forcing —»warming 
—»reduced UV-B 

Decrease in CH 4 —> less radiative forcing —* cooling 

Stratosphere 

O, formation below 18-20 km —* reduced UV-B 

0 3 formation above 18-20 km —* enhanced UV-B 

Enhanced PSC formation —» O, depletion —> enhanced UV-B 

SO x O and H,S0 4 

Troposphere 

Enhanced sulfate aerosol concentrations 

Direct radiative forcing —» cooling 

Contrail formation —> radiative forcing —» warming 

Increased cirrus cloud cover —> radiative forcing —» warming 

Modifies O, chemistry 

Stratosphere 

Modifies O, chemistry 

Soot 

Troposphere 

Direct radiative forcing —♦ wanning 

Contrail fonnation —* radiative forcing —* warming 

Increased cirrus cloud cover —*radiative forcing —* warming 

Modifies 0 3 chemistry 

Stratosphere 

Modifies O, chemistry 


In 2005, an update to the 1999 report (by Sausen 2005, based on the TRADEOFF project) 
was published to present several new values for the radiative forcing of emission components. 
Especially the RF value of contrails is strongly reduced by a factor of three to four. The fol¬ 
lowing figure shows a comparison of RF values for different components: 


26 


Aircraft RF 



Level of scientific understanding 

Good Fair Fair Fair Fair Fair Fair Poor 


Figure 2.5: Comparison of radiative forcing values for emission components (Sausen 2005) 


Out of all emission components, only contrails and nitrogen oxides (NO x ) are found to be de¬ 
pendent on the cruise altitude (Fichter 2005). CO 2 emission depends primarily on the fuel 
consumption of the aircraft. These three components will be used in the environmental impact 
model in Chapter 7. The other components have a smaller environmental impact and will not 
be considered. 


2.4.3 Contrails and Aviation-Induced Cirrus Clouds 

Contrails are the result of water vapor in the engine exhaust gasses of a jet aircraft flying at 
common cruise altitudes. The ambient temperature at this altitude is around -56° C, which 
means the water vapor freezes into tiny ice crystals (further augmented by sulfur oxides and 
soot). This formation of ice crystals forms the well-known white lines that trail after aircraft 
and can be seen from the ground. The initially thin lines can last long and spread to widths of 
more than 10 km. In the so-called airways over the North Atlantic Ocean and Europe, con¬ 
trails can cover 5% of the sky area annually. Below these airways, contrails could have a big¬ 
ger environmental impact than all greenhouse gasses combined (Whitelegg 2000). 

The presence of contrails can also induce the formation of cirrus clouds which would not nat¬ 
urally occur there. These are called aviation-induced cirrus clouds and are believed to have a 
strong warming effect on the atmosphere (Jardine 2005). In the 2005 update to the IPCC re¬ 
port, multiple studies concerning the radiative forcing of these cirrus clouds were analyzed. 
Cirrus clouds could account for the same amount of radiative forcing as all other emission 
components combined. The conclusion however is that there is too much uncertainty over the 
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actual RF value to include cirrus clouds in the total RF value of all emission components 
(Sausen 2005). This lack of scientific understanding is mentioned in Figure 2.5 as well. 

The fonnation of contrails (and thus aviation-induced cirrus clouds) depends upon the follow¬ 
ing factors: humidity, temperature, pressure, the emission index of water vapor, and the over¬ 
all propulsion efficiency of the aircraft. Furthennore, the atmosphere must be supersaturated 
with respect to ice to make ice crystal growth possible. Another significant factor is the influ¬ 
ence of latitude and the seasonal cycle, which can be seen clearly in Figure 2.6 and Figure 2.7 
(Fichter 2005). 

The TRADEOFF study following up on the aforementioned IPCC report tried to find quanti¬ 
tative results for the impact of cruise altitude changes on the global coverage of contrails and 
the radiative forcing they cause (Fichter 2005). The results of the report are summarized be¬ 
low, describing the method used and the results: 

The study applied a parameterization for line-shaped contrails by Ponater 2002. This model 
is based on the thermodynamic theory of contrail fonnation (Schmidt-Appleman theory), 
which takes into account that contrails can only form if the air is supersaturated with respect 
to ice. To calculate the actual contrail coverage from the potential contrail coverage, the dis¬ 
tance travelled was used instead of the amount of fuel used. The base case of mean contrail 
coverage was detennined by using satellite observations. And all climate change effects 
caused by contrails are measured as radiative forcing. 

Generally, a decrease in altitude results in a decrease in global contrail coverage. The rela¬ 
tionship between the two is almost linear up to a maximum decrease of 45 % in coverage at 
an altitude of 6 000 ft under the base case. On the other hand, an increase in fuel consumption 
can also be seen by flying at a lower cruise altitude (further discussed in Chapter 5 and 7). 
The report also analyzes the impact of an increase in cruise altitude. In this case there is a 
slight increase in global contrail coverage and a small decrease in fuel consumption. 
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Table 2.2: Effect of cruise altitude on contrails and fuel consumption (Fichter 2005) 


Inventory 

Flown 

Distance 

[*10 9 km/yr] 

Fuel 

Consumption 

[Tg/yr] 

Global Mean 
Contrail 
Coverage [%] 

Net RF by 
Contrails 
[mW/m 2 ] 

DLR2 

18.0 

112.2 

0.052 (distance) 
0.057 (fuel) 

2.1 (3.2) 

2.3 (3.5) 

TRADEOFF 
base case 

17.1 

111.5 

0.047 (distance) 
0.052 (fuel) 

1.9 (2.9) 

2.0 (3.1) 

TRADEOFF 

+2 kft 

17.1 

111.0 

0.050 (distance) 

2.0 (3.1) 

TRADEOFF 
-2 kft 

17.1 

114.5 

0.041 (distance) 

1.6 (2.5) 

TRADEOFF 
-4 kft 

17.1 

115.5 

0.034 (distance) 

1.3 (2.0) 

TRADEOFF 
-6 kft 

17.1 

118.0 

0.026 (distance) 

1.0 (1.6) 


The values given in the table above are all average global values. However, due to the de¬ 
pendency on latitude and the seasonal cycle, the local values of contrail coverage vary widely 
over the globe. This is visible on the following figures: 
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Figure 2.6: 


Contrail coverage variation in function of altitude changes (% in relation to the base 
case) (Fichter 2005) 
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Figure 2.7: Seasonal influence on contrail coverage (6000 ft under the base case) (Fichter 2005) 
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3 Fundamental Models 

3.1 Aerodynamics 

3.1.1 Drag Coefficient 

The drag polar expresses the relationship between lift and drag coefficient. For low-speed 
subsonic flight, a single line can represent the drag polar. When flying at higher transonic 
speeds, the drag polar consists of multiple lines. Each of these lines represent a different 
cruise Mach number. A generic high speed drag polar is given in the following figure: 



Figure 3.1 : Generic Drag polar for the transonic region (Young 2001) 

In the transonic region (where M > M crit ), the drag coefficient consists of the following 
terms: 


Cn — C D o + ACnw + K 


Dw 


r 2 

l l 


(3.1) 


With: 


K = 


1 

n ■ A ■ e 


(3.2) 


The zero lift drag coefficient C D0 depends mostly on geometric aircraft parameters and can be 
estimated with the method in Chapter 3.1.2. The wave drag increment A C Dw is highly de¬ 
pendent on the aircraft speed and is discussed in Chapter 3.1.3. The final component K . Cl is 
the lift dependent part of the drag coefficient. The Oswald factor must be known in order to 
calculate it. This is discussed in Chapter 3.1.4. 
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In the following subchapters, a lot of approximations will be made in the calculation of drag 
components. The reason being that aircraft manufactures try to keep most of their exact air¬ 
craft data classified. This means that only people in the industry can use exact data and gener¬ 
ate an accurate result. In the academic world this data is seldomly available and calculations 
must be made with approximating equations and statistical data. 


3.1.2 Zero Lift Drag Coefficient 

The zero lift drag coefficient represents the drag that is present when zero lift force is applied 
on an aircraft. It depends mostly on geometric aircraft parameters and material choice (rough¬ 
ness). It can be used to judge a design on its aerodynamic properties. 


There are a number of different options for calculating or estimating C D0 . Three different 
methods can be found in Scholz 2015. The method which is the most generally applicable is 
to calculate C D0 based on E max and the Oswald factor e (described in 3.1.4): 


Cdo — 


n ■ A ■ e 


4 -E n 


(3.3) 


With: 



(3.4) 


With: k E = 15.8 

S/Swet - 6.2 

Other methods include estimating C D0 from the wetted area of the aircraft or a more detailed 
approach by using the individual drag of aircraft components. Both methods can be found in 

Scholz 2015. 


3.1.3 Wave Drag Increment 

The wave drag increment is a term which is added to the drag coefficient equation to account 
for shock effects present at transonic speeds. Generally, the wave drag increment has a very 
small influence at low speed and only starts to gain importance when closing in on the drag 
divergence Mach number. 
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This drag divergence Mach number M dd is an important parameter in wave drag calculation. 
The so-called Boeing definition will be used here. It states that M dd corresponds to the speed 
where drag increases by 20 drag counts above the baseline condition (Blake 2009). 20 drag 
counts (CTS) equals a wave drag increment of: A C Dw — 0.0020. 

The design philosophy of Boeing and Airbus is to set M DD — M CR (Scholz 2015). But if de¬ 
sired, M dd can be calculated with a method explained in Scholz 2015. Continuing with the 
previously mentioned design philosophy, the wave drag increment can be approximated based 
on statistical data: 


AC Dw = A ■ tan (b ■ (jj—J ~ 5 ) ' cos 3 (<p 25;W ) 


(3.5) 


With: 


Merit 


B ' M dd 


tan~ 


0.002 


(3.6) 


V A ■ cos 3 (<p 25jW ) 


+ B 


An important restriction for the Mach number in Eq. (3.5) is the following: 


M < Mrri 


crit 


( 1+ ^) 


(3.7) 


A and B in these equations are constants that are based on an analysis of several aircraft. The 
values in Table 3.1 can be used when calculating one of the aircraft mentioned in the table. 
When calculating any other aircraft, the averaged values in Table 3.2 should be used. The Ex¬ 
cel tool uses the constants from Table 3.2 as default. Table 3.1 clearly shows that the parame¬ 
ters can vary widely between different aircraft models, which means that the average statisti¬ 
cal values used in the Excel tool make the model slightly less accurate. 

Table 3.1: Statistically derived values for A and B for a number of aircraft (Scholz 2015) 



A 320-200 

B727-200 

B737-800 

C-130H 

BAe 146 

A 

0.000885 

0.000766 

0.001171 

0.001201 

0.001765 

B 

3.734 

5.257 

3.543 

3.126 

3.457 


Table 3.2: Statistically derived values for A and B for any aircraft type (Scholz 2015) 



Result 

A 

0.001272 

B 

3.477 
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3.1.4 Oswald Factor 

The Oswald factor or Oswald efficiency number is a value which, together with the aspect ra¬ 
tio, defines the rate of change of the drag coefficient with a varying lift coefficient. It is used 
to calculate the induced drag. 

For an estimation of the Oswald factor e a method from Scholz 2015 will be used. Because 
the estimated e is used in the calculation of C D0 (see Chapter 3.1.2), the method to estimate e 
should not rely on the value of C D0 . Thus, ‘method 1’ is used without the input of C D0 . The 
Oswald factor is composed of the following factors: 

£ ®theo ' k e ,F ' ke,D0 ' ^e,M (3-8) 

The limiting factor here is the correction factor for losses due to the compressibility effects on 
induced drag: k e M . It is dependent on the flight Mach number and is calculated as follows: 

/ M \ be 

k-e,M — ®-e ' I 77 1 ] T C e (3.9) 

V v *comp ) 


With: 


b e = 10.82 
c e = 1 
Mcomp — 0.3 


And: 


a e 


-1 


M + 0.08 

. M C omp 



Equations for the other factors can be found in Scholz 2015. 


(3.10) 


3.2 Thrust Specific Fuel Consumption 

The TSFC is defined as the mass of fuel burned per unit of time, divided by the thrust. It is an 
important perfonnance parameter for a jet engine. A smaller TSFC means a better effeciency 
of the engine (Young 2018). 

In theory TSFC is dependent on altitude (overall engine efficiency), thrust and speed. The 
overall engine efficiency depends partially on the difference between maximum internal air 
temperature and outside air temperature. Implying that overall engine efficiency increases 
with rising altitude. However, most commercial jet aircraft cruise in the stratosphere where 
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the temperature is invariable with altitude (see Chapter 3.4). In the altitude range of cruise 
flight for these aircraft, the variation in TSFC varies between 1% and 2%. Thus, the 
dependency on altitude can be neglected according to Bensel 2018. 

The dependency on thrust can also be considered small enough to neglect. This is proven by 
Bensel 2018 by retrofitting data from 4 different turbojet engines to visualize a general trend 
of very small variation in TSFC with varying amounts of thrust during cruise (see Figure 3.2). 
With thrust variatons of up to 25 %, the change in TSFC is less than 5 %. This conclusion is 
shared by Young 2018. 

As a result of this, TSFC would only be dependent of aircraft speed. An important note to 
make here is that in the study by Bensel the aircraft are considered to be flying at common 
cruise altitudes for commercial jet aircraft. For the scope of this thesis this is not always the 
case. This means that the altitude-dependency will be taken into account in the TSFC model 
in Chapter 5. 
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Figure 3.2: Design TSFC in function of design thrust (Bensel 2018) 
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The TSFC variation with speed can be modelled linearly according to the model from Roux 
2002 and Scholz 2017 or the non-linear model by Herrmann 2010. Taking into consideration 
that most modem turbofan engines have a relatively high bypass ratio, Bensel 2018 concluded 
that the non-linear Herrmann model fits closer to the actual TSFC characteristics of a high- 
BPR turbofan engine. The input variables are also easier to obtain than those for the linear 
model. The non-linear model does include the impact of altitude on the TSFC as well. 
Therefore, the Herrmann model will be used for all further TSFC calculations in this thesis. 

All equations used in the model by Herrmann 2010 can be found in Appendix A. 

The input variables are: 

- Cruise Mach number 

- Cruise altitude 

- Bypass ratio 

- Take-off thrust (one engine) 

- Overall pressure ratio (if known) 

- Turbine entry temperature 

As an output TSFC in kg/Ns is generated. 


3.3 Direct Operating Costs 

The DOC is calculated with a method from Scholz 2015 based on the AEA1989 method. Es¬ 
sentially, direct operating costs only contain aircraft-related costs and can be divided into 7 
categories and several subcategories. These categories are then added up to fonn one DOC 
value, expressed as an annual DOC cost for one aircraft. However, the more accurate way to 
compare multiple aircraft with different design and flight missions is by using the seat-mile 
cost: 
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° s,m 


c, 


a/c,a 
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DOC 


(3.11) 
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The full list of categories and subcategories: 

• Depreciation C DEP 

• Interest C INT 

• Insurance C INS 

• Fuel C F 

• Maintenance C M : 

o Airframe maintenance C MAF 
o Power plant maintenance C MPP 

• Crew C c : 

o Cockpit crew C c co 
o Cabin crew C c ca 

• Fees and charges C FEE : 

o Landing fees C FEE LD 
o ATC or navigation fees C FEENAV 
o Ground handling fees C FEE GND 

All equations and necessary data can be found in Scholz 2015 or in the Excel tool as an ap¬ 
plied model. 


3.4 Atmosphere 

From a flight mechanics point of view the atmosphere is often modelled by the standard at¬ 
mosphere. This is an idealized, steady-state representation of the Earth’s atmosphere con¬ 
structed by different organizations: ISO with ISO 2533:1975, ICAO with the International 
Standard Atmosphere and the United States government with the US Standard Atmosphere. 
Up to an altitude of 32 km all three models are identical (NASA, 1976). Within the scope of 
this thesis, the layers of interest are the two lowest layers of the Earth’s atmosphere: the trop¬ 
osphere and stratosphere, which are divided by the tropopause. 

Note that the further mentioned altitude of the tropopause of 11 000 m (36 089 ft) is an aver¬ 
age value used in the standard atmosphere model. The actual altitude of the boundary between 
troposphere and stratosphere varies significantly between the Equator and the poles. 


The following mathematical approach is summarized from Young 2001: 
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1) Troposphere: 

This is the region which spans from sea level up to the tropopause at 11 000 m (36 089 ft). 
The temperature is assumed to be exactly 288.15 K at sea level and to decrease linearly with 
0.0065 K/m for increasing altitude. 


Temperature, pressure, and density are generally given as ratios of the ambient condition to 
the condition at sea level (denoted with subscript 0). 
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6 = — with-. T — T 0 — L ■ H (3.12) 
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(3.13) 


(3.14) 


With: T 0 =288.15 K 

L = 0.0065 K/m 
Po = 101 325 Pa 
Po = 1.2250 kg/m 3 
g = 9.80665 m/s 2 
R = 287.053 m 2 /s 2 /K 

2) Stratosphere: 

This is the region above the troposphere. The focus here is the part of the stratosphere from 
the tropopause at 11 000 m up to 20 000 m (65 617 ft). Within this zone it is assumed that the 
temperature is a constant 216.65 K. 

In this case temperature, pressure, and density are given as ratios of the ambient condition to 
the conditions at tropopause level (denoted with subscript t). This simplifies into two equa¬ 
tions: 


T = T t — 216.65 K (3.15) 

— — — = ^ H ~ H T> 

Pt Pt 


(3.16) 
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With: p t = 22 657 Pa 

p t = 0.36392 kg/m 3 
H t = 11 000 m 
^ = 9.80665 m/s 2 
= 287.053 m 2 /s 2 /K 


The following figure illustrates the behavior of the three main atmospheric parameters in 
function of altitude up to 15 000 m: 


Temperature Pressure 




Density 



Density (kg/m3) 


Figure 3.3: Distribution of atmospheric parameters in function of altitude 


3.5 Equivalent CO 2 Mass 


The equivalent CO 2 mass or m C02ieq is a metric used to express the environmental impact of 
engine emissions in a quantitative way. Introduced in Schwarz 2011, it combines the effects 
of the three biggest emission contributors: CO 2 , NO x , and aircraft-induced cloudiness (AIC). 
AIC includes contrails and aircraft-induced cirrus clouds (see Chapter 2.4.3). The equation is 
made as generic as possible by dividing equivalent CO 2 mass by the number of seats on the 
aircraft. This results in units of kg CO 2 per NM flown, per seat: 
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(3.17) 
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To get a better understanding of this equation, it is displayed with its units only: 

kg C02 ___ 

NM ■ n sea t 

kg C02/kg fuel ■ kg fuel/NM i kg NOx/kg fuel ■ kg fuel/NM kg C02 NM kg C02 

-'1H-' --1-- (o.loj 

Tiseat Yiseat kg NOx NM w Yi sea t NM 

The equation is made up of two types of u nkn own parameters: Emission indices (El) and 
characterization factors (CF). The emission index is the mass amount of the specie emitted per 
mass amount of fuel burned. The emission index of CO 2 is equal to 3.16 kg CO 2 per kg of 
fuel burned and is independent of the aircraft’s altitude (Schwarz 2011). The emission index 
of NOx is a variable value that does depend on the altitude. It is calculated using the Boeing 
fuel flow method 2 (Baughcum 1996). 

The Boeing fuel flow method 2 (BFFM2) is a method for calculating the emission index of 
NOx based on fuel flow and emission index data from the ICAO Aircraft Engine Emissions 
Databank (ICAO 2019). The method is adjusted for atmospheric effects at different altitudes 
and represents a relationship between fuel flow and the emission index of NO x . Essentially, it 
starts from the actual fuel flow of the considered aircraft (calculated from its fuel consump¬ 
tion and TAS) and corrects it for the reference condition. Then, El of NO x and fuel flow val¬ 
ues from the ICAO databank are corrected by a factor defined by Boeing and plotted in a log- 
log plot for several different flight stages. The next step is to derive an equation for El NO x in 
function of fuel flow by curve fitting the log-log plot. The final step is to correct the El NO x 
value for actual atmospheric conditions. The method is described in more detail in Appendix 
B and in an applied fonn for the Excel tool in Chapter 4.7. 

The characterization factor (CF) may be considered as a conversion factor from NO x or AIC 
emissions to the equivalent CO 2 emission. Both are altitude dependent by introducing a forc¬ 
ing factor s. The result has units of kg of CO 2 per kg of NO x and kg of CO 2 per NM flown re¬ 
spectively: 


C ^'midpoint, N Ox 


SGTP, 


03s,100 


SGTP C0 2 ioo 


S 03s(H ) + 


SGTP, 


03L,100 


SGTP C0 2 100 


s 03l(H) + 


SGTP 


CH 4,100 


SGTP C0 2,100 


s CH4 (H) (3.19) 


r _ SGTP cont raiis, 100 SGTP clrrus 100 

c r midpoint,AIC <:rTP b contrails\ n ) ' crTP ^cirrus\ n J 

2,100 2,100 


The SGTP factor stands for sustained global temperature potential. Environmental impact is 
generally measured in global wanning potential (GWP), which uses radiative forcing to ex¬ 
press the impact of a system. Radiative forcing however is not the best metric for estimating 
the impact of a single flight (see Chapter 2.4.2). Instead SGTP is used as an alternative metric. 
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It expresses an average temperature increase or decrease (over 100 years) in Kelvin per unit of 
the respective specie. Units are K per relevant unit for each of the emission contributors (per 
kg CO 2 , per kg NO x or per NM for AIC). SGTP values can be found in the following table: 


Table 3.3: SGTP values (Schwarz 2009) 


Species 

SGTPi.ioo 

CO 2 (K/kg CO 2 ) 

3,58 • 

10-14 

Short O 3 (K/kg NOx) 

7,79- 

10- 12 

Long O 3 (K/NO x ) 

-9,14 

• 10- 13 

CH 4 (K/kg NOx) 

-3,90 

■ 10- 12 

Contrails (K/NM) 

2,54 • 

10- 13 

Contrails (K/km) 

1,37- 

10- 13 

Cirrus (K/NM) 

7,63 ■ 

10- 13 

Cirrus (K/km) 

4,12 • 

10- 13 


The forcing factor 5 introduces the altitude dependency in Eq. (3.19) and (3.20). The forcing 
factors for short ozone, long ozone and methane increase with increasing altitude. The AIC 
forcing factor reveals a clear peak around the common cruise altitude for jet aircraft and turns 
to almost zero at lower altitudes (see Figure 3.4). Exact values for all forcing factors are ex¬ 
tracted in Van Endert 2017 and can be found in the Excel tool and in Appendix C. 



Figure 3.4: Forcing factor s in function of altitude (Schwarz 2011) 
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4 Excel Tool 

4.1 General 

In the four following chapters calculations will be made to come up with results for the differ¬ 
ent models. To keep an overview, all this is done in an Excel tool which will be added to this 
thesis. This chapter delivers a brief overview of the tool. Each subchapter, starting from Chap¬ 
ter 4.4, is dedicated to one sheet of the Excel tool. Most of the calculations in the tool are 
done by macros in the form of buttons. All these macros can be executed separately, there is 
no need to calculate them in any order. For ease of use, the following color code is used con¬ 
sistently throughout the Excel tool: 


Table 4.1 : Color code used in the Excel tool 


User input value 

Bold black on green background 

Calculated input value 

Green on white background 

Important result 

Black on orange background 

Result represented in plot 

Black on blue background 


An important note for the user is that the tool was written in an Excel version that uses a 
comma as a decimal separator instead of a dot. This means that all results acquired from the 
tool will be displayed with a comma as decimal separator. 


4.2 Case Study: Airbus A320-200 

As a representative example the Airbus A320-200 (weight variant 000) will be used in all 
models. This is a popular aircraft at the time of publication an has been for some time, mean¬ 
ing that a great amount of aircraft data can be found relatively easy. Furthermore, it is a single 
aisle aircraft. This is the group of aircraft with the biggest global market share by far right 
now and probably will be in the near future (69 % of all passenger aircraft and growing) ac¬ 
cording to Boeing 2018. All aircraft data used in the models is found in Airbus 2005a and 
Jackson 2007. 

4.3 Macro for Calculation 

All calculations are done with macros programmed in Visual Basic in Excel. This is a quick 
example of the standard loop code that is used to calculate the matrix style output values. In 
this way the user can insert additional restrictions or alter the code. The if statement inside the 
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loop is meant as a safety to prevent output values being displayed for Mach numbers greater 
than the maximum operating Mach number of the aircraft. 


Sub Calculate_fuelcons() 

For i = 52 To 71 

j = 3 To 13 

Sheets("Inputs_Outputs") .Cells (i, 2) .Copy 
Sheets("Inputs_Outputs") .Cells (3, 10) 

Sheets("Inputs_Outputs") .Cells (51, j) .Copy 
Sheets("Inputs_Outputs") .Cells (2, 10) 

If Sheets("Inputs_Outputs").Cells(51, j).Value < 

Sheets("Inputs_Outputs") .Cells (14, 2).Value Then 

Sheets("Fuel") .Cells (39, 9) .Copy 

Sheets("Inputs_Outputs") .Cells (i, j) .PasteSpecial xlPasteValues 
Else 

Sheets("Inputs_Outputs").Cells(i, j).Value = "" 

End If 
j 

Next i 

End Sub 


4.4 Inputs and Outputs 

The inputs are divided into categories: aircraft parameters, engine parameters, flight parame¬ 
ters, and constants. 

- Aircraft parameters: these parameters are used for all outputs and should be filled in as 
thorough as possible. 

- Engine parameters: are used for TSFC calculation and for calculating maintenance 
costs and engine price in the DOC. 



43 


- Flight parameters: contain the two main inputs in this model: Mach number and alti¬ 
tude. All other parameters (except range) change with varying Mach number and alti¬ 
tude. The range is needed for DOC calculation. 

- Constants: the values are taken from Young 2001 and can be left as is. 

Below the input categories, all outputs are shown. From top to bottom are four different cate¬ 
gories and from left to right in general the main parameter followed by the most influential 
parameters for each category and/or the nonnalized value of the main parameter. The follow¬ 
ing lay-out is used: 


Table 4.2: Lay-out of the outputs in the Excel tool 


Lift coefficient 

Drag coefficient 

Oswald factor 

Wave drag incre¬ 
ment 

Fuel consumption 

TSFC 

Aerodynamic effi¬ 
ciency 

Fuel consumption 
normalized 

Seat-mile cost 

Annual DOC per air¬ 
craft 

Flight time 

Seat-mile cost 
normalized 

Environmental im¬ 
pact 

Equivalent CO 2 mass 

EINOx 

Equivalent CO 2 
mass normalized 

Combined values 





Note that the input Mach number goes up to 0.9, but the displayed results are limited by the 
chosen aircraft’s maximum operating Mach number. All outputs for a higher Mach number 
will be left blank. 

At the bottom of the sheet, four weighting factors can be found. Both the environmental im¬ 
pact model and the combined model are made up of different values with a weighting factor 
each. The environmental impact model is composed of resource depletion (fuel consumption) 
and equivalent CO 2 mass. The combined model is composed of seat-mile cost and environ¬ 
mental impact. All weighting factors are set to 0.5 by default, which means that each compo¬ 
nent has the same importance. The user can change the influence of each component on the 
total result by altering the weighting factors. 

All outputs that are multiplied with a weighting factor are normalized values. Calculation re¬ 
sults are first normalized between 0 and 1, then multiplied with their respective weighting fac¬ 
tor and finally added up to get a combined value between 0 and 1: 

combined value = 


(wf 1 ■ normalized val. 1) + (wf 2 ■ normalized val. 2) + ••• 


(4.1) 




44 


4.5 Fuel 

The fuel consumption sheet is used to calculate all aerodynamic parameters on the left: Os¬ 
wald factor, wave drag increment, lift- and drag coefficient, and the aerodynamic efficiency. 
On the right side TSFC is calculated with the method from Herrmann 2010. 

In this sheet no user input is necessary. Its main purpose is to elaborate the calculation meth¬ 
ods for each parameter of the fuel consumption calculation. 


4.6 DOC 

This sheet is used to calculate the annual DOC and the seat-mile cost. It also generates a visu¬ 
alization of the weight of each category from the total DOC by using a pie chart. This is locat¬ 
ed under the total annual DOC value. 

At the bottom of the sheet is a table with matching bar graph. Here the user can enter four dif¬ 
ferent flight scenarios (Mach and altitude) and compare the DOC distribution for each scenar¬ 
io. 


The following user input is required for this sheet: 

1. The current year of calculation must be entered for a correct inflation factor. 

2. A selection must be made between a short/medium-range aircraft or a long-range air¬ 
craft. 

3. The current price for jet A-l fuel (kerosene) must be entered from EIA 2019 or an¬ 
other source (must be FOB price). 


4.7 Environmental 

In this sheet the environmental impact of the aircraft is calculated. This is the combination of 
equivalent CO 2 mass per Nm flown per seat and resource depletion (as fuel consumption). At 
the bottom of the sheet the user has the option to compare four flight scenarios (Mach number 
and altitude) and the corresponding impact of each emission component. 

There is some user input required in this sheet to calculate the emission index for NOx. This 
is done by using the Boeing fuel flow method 2 (BFFM2) (Baughcum 1996). Based on the 
fuel flow and EI_NO x values from the ICAO Aircraft Engines Emission Database (ICAO 
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2019). The full method is introduced in Chapter 3.5 and fully explained in Appendix B. In the 
Excel tool the following steps must be taken by the user: 


1. The unadapted fuel flow and EI_NO x values for T/O, C/O, App and Idle must be ob¬ 
tained from ICAO 2019 for the appropriate engine and filled in in the following table 
(green columns only): 


Table 4.3: ICAO engine data (ICAO 2019) 


Flight stage 

Wf_unadapted 

r 

Wf_adapte EI_NOx 


T/O 

1,166 

1,01 

1,17766 

28,7 

C/O 

0,961 

1,013 

0,973493 


23,3 

App 

0,326 

1,02 

0,33252 


10 

Idle 

0,107 

1,1 

0,1177 

4,3 


2. The graph accompanying the table shows the correlation between adapted fuel flow 
and the corresponding EI_NO x values, including a trendline for this correlation (see 
Figure 4.1). The trendline function should be filled in in cell C45 (also changing the x 
variable in the function to cell C43) to get an uncorrected EI_NO x value for the cruise 
stage. With some further calculations following the Boeing fuel flow method 2, the final 
corrected EI_NO x value is presented in cell C51. 


100 

y = 22,412x +1,9994 

R 2 = 0,998 

. 

• 

*. 



0,1 1 10 

log Wfadapted 


Figure 4.1: Correlation between adapted fuel flow and EI_NO x 


4.8 Flight Time 

The flight time sheet is used to generate an approximating model for a generic passenger jet 
aircraft. No user input is required. BADA data (Eurocontrol 1998) is presented for four air¬ 
craft: Airbus A320-200, Airbus A340, Boeing 737-200, and Boeing 747-400. This data is pre¬ 
sented in graphs and is used to calculate an averaged model. The averaged model is then 
compared to the results generated by using Eq. 6.1 (a simple flight time calculation). All cor¬ 
responding graphs can be found in Appendix D as well. 
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4.9 Extra Information 


This sheet includes a linear interpolation model for the forcing factor data (see Appendix C), 
the Airbus A320-200 payload-range diagram with a simple fuel consumption calculation (see 
Chapter 5.1), and the Airbus A320-200 drag polar based on data from the fuel consumption 
sheet. No user input is required in this sheet. 
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5 Flight Parameter Optimization: Minimum Fuel 
Consumption 

In the following four chapters all optimization results are generated for an Airbus A320-200 
(weight variant 000) as a representative example but the models are made for any passenger 
jet aircraft. 


5.1 General 

In this chapter a simple model will be designed and utilized to calculate an aircraft’s fuel con¬ 
sumption. It is mostly based on aerodynamic parameters and TSFC, in function of the air¬ 
craft’s cruise Mach number and altitude. 

As mentioned in Chapter 2.1, minimum fuel consumption during cruise is related with maxi¬ 
mum specific air range r a . The equation for the SAR (Eq. (2.6)) can be altered so it contains 
the Mach number: 


V ■ E M-a-E 

r a — -=- ( 5 . 1 ) 

c ■ m - g c ■ m - g 

The SAR is expressed in units of km/kg fuel, which means the inverse f = l/r a can be de¬ 
fined as the fuel consumption in units of kg fuel per km (or NM) flown. Nautical mile will be 
used throughout the model instead of kilometer. The mass of fuel is used instead of the vol¬ 
ume because the density of fuel changes with altitude. Fuel consumption can be calculated 
with the following equation: 


c -m - g 
M-a-E 


(5.2) 


With g as a constant, m depending on the aircraft and mission, and all other parameters de¬ 
pending on flight parameters. However, to make the model as generic as possible / will be 
divided by m (which is the MTOW in this case). This results in the following equation in 
units of kg fuel per NM flown, per kg of MTOW: 


/mto 


w 


c-g 

M-a-E 


(5.3) 


An alternative way to Eq. (5.2) to calculate fuel consumption is by calculating the specific air 
range with the aircraft’s payload-range diagram. This should be used as an estimation rather 
than a real calculation. Or it could be used to check the order of magnitude of a calculated fuel 
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consumption result. The drawback of this option is that the required values for range and pay- 
load mass are read from a diagram, which comes with a certain level of inaccuracy. 


The method is based on calculating specific air range by using maximum range and range for 
maximum payload with their respective payloads. The specific air range is defined as: 


R z R± 
m 1 — m 2 


(5.4) 


With R and m values depicted on the following example of a payload-range diagram: 



Figure 5.1: Illustrative example of a payload-range diagram (Van Endert2017) 

With: Ri — range at maximum payload 

R 2 — maximum range 
m 1 — maximum payload 
m 2 — payload at maximum range 

An aircraft’s payload-range diagram is usually found in the publicly available “documents for 
airport planning”. As an example, the Airbus A320 document for airport planning (Airbus 
2005a) can be used for the Airbus A320-200 payload-range diagram. 

Continuing with Eq. (5.2) to calculate fuel consumption: M and C L were expected to be the 
inputs in this model according to the task sheet. These two parameters are mutually depend¬ 
ent, with the altitude H as the link between them. The lift coefficient C L is defined as: 
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C L 


2 ■ m ■ g 
p ■ M 2 ■ a 2 ■ S 


(5.5) 


This equation implicates that the relationship between M and C L depends on p, a and S. The 
first two parameters are in turn dependent on the altitude H. The dependency of the density p 
is given in Chapter 3.4 and for a in Chapter 5.2.1. Using the following set of equations: 


a = yfy ■ R-T 


T — T 0 — L ■ H 


P = Po ■ 


L ■ H 


T 0 J 


RL 


(5.6) 

(3.12) 

(3.14) 


Eq. (5.5) can be solved for the altitude H: 


To_ f{ 2-m-g \ R ' L/g 

L \Vp 0 ' T 0 ■ M 2 ■ y ■ R ■ S ■ C L ) 


(5.7) 


The problem with this model is that the mutual dependency of M and C L is also influenced by 
the wing surface area S. And because S is different for every aircraft, the model would be¬ 
come less generally applicable. For this reason, the choice has been made to deviate from the 
original task description and use the Mach number M and the altitude H as inputs in all mod¬ 
els. 


5.2 Important Parameters 

The most important parameters for calculating fuel consumption according to Eq. (5.3) are 
displayed in the following diagram. Inputs are green, important intermediate parameters blue 
and the final result in grey. 
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Figure 5.2: Flowchart of fuel consumption calculation 


5.2.1 Speed of Sound 

The speed of sound a is dependent on the ambient temperature and can be calculated as fol¬ 
lows: 


a — yjy ■ R - T 


(5.8) 


With: 7 = 1.4 

R = 287.053 m 2 /s 2 K 
T = ambient temperature in K 

Due to the change of temperature with varying height (see Chapter 3.4), the speed of sound is 
dependent on cruise altitude. 


5.2.2 Aerodynamic Efficiency 

The aerodynamic efficiency E is the ratio between lift and drag, and can be written in several 
different ways: 


L _ 2 ■ p ■ S ■ V 2 ■ C L _ C L 
D ~ 2 ■ p ■ S ■ V 2 ■ C D ~ Cb 


(5.9) 
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By introducing Eq. (3.1), E can be written as: 


E = 


C, 


Cdo 3“ AC, 


Dw 


+ K-C} 


(5.10) 


With: 


K 


n ■ A ■ e 


(3.2) 


A model must be used to calculate C D0 , A C Dw , and e. This is accomplished by using different 
methods from Scholz 2015 (see Chapter 3.1). The lift coefficient C L is defined as: 


C L 


2 ■ m ■ g 
p ■ M 2 ■ a 2 ■ S 


(5.5) 


This equation and Eq. (3.2) imply that the aerodynamic efficiency E is partially dependent on 
geometric aircraft design parameters including the aspect ratio A and the surface area S. 


5.2.3 Mach Number 

The Mach number M is considered as an input in this model. It can be derived from true air¬ 
speed and local speed of sound if need be: 


V 

M — — (5.11) 

a 


5.2.4 Specific Fuel Consumption 

The calculation of the specific fuel consumption c is accomplished with the non-linear meth¬ 
od provided in Herrmann 2010. This method is introduced and further elaborated in Chapter 
3.2. All necessary equations can be found in Appendix A. 


5.3 Results Fuel Consumption 

The results are presented as generic as possible. By calculating fuel consumption as kg of fuel 
per NM flown and per kg of aircraft MTOW, aircraft of different sizes and for different flight 
missions can be compared. Even an aircraft to be designed can be used without knowledge of 
its MTOW. On the other hand, the most common units used in relevant literature are kg of 
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fuel per 100 km flown per seat (comparable to the units used in the automotive industry). This 
representation will be included as well. 


5.3.1 TSFC 


Thrust specific fuel consumption is represented as kg of fuel used per Newton of thrust per 
second. And this in function of Mach number and altitude. The lower the value is, the better. 


Mach number 



0,4 

0,45 

0,5 

0,55 

0,6 

0,65 

0,7 

0,75 

0,8 

3000 

l,08E-05 

l,15E-05 

1.22E-05 

l,3E-05 

l,39E-05 

l,49E-05 

1.6E-05 

l,73E-05 

l,87E-05 

3500 

l,08E-05 

l,15E-05 

l,22E-05 

l,3E-05 

l,39E-05 

l,48E-05 

l,59E-05 

l,71E-05 

l,85E-05 

4000 

l,08E-05 

l,15E-05 

l,22E-05 

l,3E-05 

l,38E-05 

l,47E-05 

l,58E-05 

l,7E-05 

l,83E-05 

4500 

l,08E-05 

l,15E-05 

1.22E-05 

l,29E-05 

l,37E-05 

l,47E-05 

1.57E-05 

l,68E-05 

l,81E-05 

5000 

l,08E-05 

l,15E-05 

l,21E-05 

l,29E-05 

l,37E-05 

l,46E-05 

l,56E-05 

l,67E-05 

l,79E-05 

5500 

l,08E-05 

l,15E-05 

l,21E-05 

l,28E-05 

l,36E-05 

l,45E-05 

l,55E-05 

1.65E-05 

l,78E-05 

6000 

l,08E-05 

l,14E-05 

1.21E-05 

l,28E-05 

l,36E-05 

l,44E-05 

l,54E-05 

l,64E-05 

l,76E-05 

6500 

l,08E-05 

l,14E-05 

l,21E-05 

l,28E-05 

l,35E-05 

l,44E-05 

l,53E-05 

l,63E-05 

l,74E-05 

7000 

l,08E-05 

l,14E-05 

l,21E-05 

l,28E-05 

l,35E-05 

l,43E-05 

l,52E-05 

l,62E-05 

l,73E-05 

7500 

l,08E-05 

l,14E-05 

1.21E-05 

l,27E-05 

l,35E-05 

l,42E-05 

l,51E-05 

l,61E-05 

l,71E-05 

8000 

l,08E-05 

l,14E-05 

l,20E-05 

l,27E-05 

l,34E-05 

l,42E-05 

l,50E-05 

1.60E-05 

l,70E-05 

8500 

l,08E-05 

l,14E-05 

l,20E-05 

l,27E-05 

l,34E-05 

l,41E-05 

l,50E-05 

l,59E-05 

l,69E-05 

9000 

l,08E-05 

l,14E-05 

1.20E-05 

l,26E-05 

l,33E-05 

l,41E-05 

l,49E-05 

l,58E-05 

l,67E-05 

9500 

l,09E-05 

l,14E-05 

l,20E-05 

l,26E-05 

l,33E-05 

l,40E-05 

l,48E-05 

l,57E-05 

l,66E-05 

10000 

l,09E-05 

l,14E-05 

l,20E-05 

l,26E-05 

l,33E-05 

1.40E-05 

l,48E-05 

l,56E-05 

l,65E-05 

10500 

l,09E-05 

l,14E-05 

l,20E-05 

l,26E-05 

l,32E-05 

l,39E-05 

l,47E-05 

l,55E-05 

l,64E-05 

11000 

l,09E-05 

l,14E-05 

l,20E-05 

l,26E-05 

l,32E-05 

l,39E-05 

l,46E-05 

l,54E-05 

1.63E-05 

11500 

l,09E-05 

l,14E-05 

l,20E-05 

l,26E-05 

l,32E-05 

l,39E-05 

l,46E-05 

l,54E-05 

l,63E-05 

12000 

l,09E-05 

l,14E-05 

l,20E-05 

l,26E-05 

l,32E-05 

l,39E-05 

l,46E-05 

l,54E-05 

l,63E-05 

12500 

l,09E-05 

l,14E-05 

l,20E-05 

l,26E-05 

l,32E-05 

l,39E-05 

l,46E-05 

l,54E-05 

l,63E-05 


Figure 5.3: Results TSFC in kg of fuel per N of thrust per second 


TSFC increases with increasing Mach number and decreasing altitude (except for the case of 
Mach 0.4, which stays almost constant). TSFC depends on the overall engine efficiency rj 0 
which is a combination of five engine component efficiencies according to the model by 
Herrmann 2010 (see Appendix A): inlet, fan, compressor, turbine, and nozzle. Four of these 
(all except inlet efficiency) are directly influenced by the cruise Mach number in a negative 
way. The gas generator efficiency equation and the final equation for calculating the TSFC 
are also directly influenced by the Mach number. Altitude influences the TSFC model by 
means of the ambient temperature. Ambient temperature has an influence on the gas generator 
function and on the final TSFC equation. 


Apart from ambient temperature, altitude influences the TSFC in another more profound way 
which is possibly not as visible by representing results in the style of Figure 5.3. The speed of 
sound varies with altitude and by using Mach number instead of TAS as columns, the TAS 
does not remain constant throughout one column. An increase in altitude means a decrease in 
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TAS for a certain Mach number. This explains the trend of lower TSFC at higher altitude: the 
TSFC is lower because the TAS is lower for an equal Mach number. 

The relationship between TSFC and TAS is displayed in Figure 5.4. Starting from the bottom 
left corner, the three different shapes which are on almost the same horizontal level represent 
the TSFC at one Mach number. Moving to a higher speed, the TSFC increases with decreas¬ 
ing altitude for an equal Mach number. At a TAS of 0 knot, the TSFC is not equal to zero 
(engine on idle). The lower the TAS goes to 0, the closer together all TSFC values are, inde¬ 
pendent of altitude. This effect is visible in the results on Figure 5.3 as well. 
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Figure 5.4: TSFC in function of true airspeed for three different altitudes 


The calculated output for the Airbus A320 with CFM56-5B4 engines at regular flight condi¬ 
tions (Mach 0.78 at an altitude of 11 500 m) is 1.59 ■ 10~ 5 kg/Ns, which corresponds with the 
value range found in Jenkinson 1999 for the CFM56 engine family: 1.54 ■ 10~ 5 — 1.89 ■ 
10~ 5 kg/Ns. This validates the accuracy of the TSFC model to a certain extent. 
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5.3.2 Aerodynamic Efficiency 


The aerodynamic efficiency E is a dimensionless parameter. The higher it is, the better. 



Mach number 


0,4 

0,45 

0,5 

0,55 

0,6 

0,65 

0,7 

0,75 

0,8 


3000 

18,59 

19,47 

19,46 

18,89 

17,96 

16,71 

15,38 

13,90 

11,88 


3500 

18,19 

19,32 

19,55 

19,16 

18,39 

17,23 

15,95 

14,49 

12,45 


4000 

17,73 

19,08 

19,55 

19,37 

18,75 

17,71 

16,50 

15,07 

13,02 


4500 

17,22 

18,77 

19,47 

19,51 

19,06 

18,14 

17,02 

15,64 

13,59 


5000 

16,65 

18,39 

19,31 

19,56 

19,29 

18,52 

17,50 

16,18 

14,16 


5500 

16,05 

17,94 

19,06 

19,52 

19,44 

18,84 

17,93 

16,70 

14,71 


6000 

15,41 

17,42 

18,73 

19,39 

19,51 

19,08 

18,31 

17,18 

15,25 


6500 

14,75 

16,85 

18,32 

19,18 

19,49 

19,24 

18,62 

17,60 

15,76 


7000 

14,07 

16,24 

17,84 

18,87 

19,37 

19,30 

18,85 

17,96 

16,22 

a; 

-Q 

7500 

13,38 

15,58 

17,30 

18,48 

19,16 

19,28 

18,99 

18,25 

16,64 

13 

8000 

12,68 

14,90 

16,69 

18,01 

18,85 

19,15 

19,03 

18,46 

17,00 

■*-> 

< 

8500 

11,98 

14,19 

16,04 

17,47 

18,45 

18,92 

18,98 

18,58 

17,29 


9000 

11,29 

13,47 

15,35 

16,86 

17,97 

18,60 

18,83 

18,60 

17,49 


9500 

10,61 

12,74 

14,63 

16,20 

17,41 

18,18 

18,57 

18,51 

17,60 


10000 

9,95 

12,01 

13,89 

15,50 

16,79 

17,68 

18,22 

18,32 

17,61 


10500 

9,30 

11,29 

13,13 

14,76 

16,10 

17,10 

17,77 

18,03 

17,52 


11000 

9,38 

11,37 

13,22 

14,85 

16,19 

17,17 

17,82 

18,07 

17,54 


11500 

8,75 

10,66 

12,46 

14,08 

15,46 

16,53 

17,29 

17,68 

17,35 


12000 

8,15 

9,96 

11,71 

13,31 

14,70 

15,83 

16,69 

17,21 

17,07 


12500 

7,58 

9,30 

10,98 

12,54 

13,94 

15,10 

16,03 

16,65 

16,69 


Figure 5.5: Results aerodynamic efficiency 


The aerodynamic efficiency E has an optimum zone that follows a rising Mach number - ris¬ 
ing altitude pattern. E is defined as E — C L /C D and this ratio clarifies the pattern. Looking at 
the following two figures for lift coefficient and drag coefficient, the optimum zone for E is 
associated with moderate values for both coefficients (yellow zone). This is because a high lift 
coefficient has a high drag coefficient as a result and vice versa, which alters the ratio be¬ 
tween the two. 
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Mach number 



0,4 

0,45 

0,5 

0,55 

0,6 

0,65 

0,7 

0,75 

0,8 

3000 

0,818 

0,647 

0,524 

0,433 

0,364 

0,310 

0,267 

0,233 

0,205 

3500 

0,873 

0,689 

0,558 

0,461 

0,388 

0,330 

0,285 

0,248 

0,218 

4000 

0,931 

0,736 

0,596 

0,492 

0,414 

0,353 

0,304 

0,265 

0,233 

4500 

0,994 

0,785 

0,636 

0,526 

0,442 

0,376 

0,325 

0,283 

0,248 

5000 

1,062 

0,839 

0,680 

0,562 

0,472 

0,402 

0,347 

0,302 

0,266 

5500 

1,136 

0,898 

0,727 

0,601 

0,505 

0,430 

0,371 

0,323 

0,284 

6000 

1,216 

0,961 

0,778 

0,643 

0,541 

0,461 

0,397 

0,346 

0,304 

6500 

1,303 

1,030 

0,834 

0,689 

0,579 

0,493 

0,425 

0,371 

0,326 

7000 

1,397 

1,104 

0,894 

0,739 

0,621 

0,529 

0,456 

0,397 

0,349 

7500 

1,500 

1,185 

0,960 

0,793 

0,667 

0,568 

0,490 

0,427 

0,375 

8000 

1,612 

1,274 

1,032 

0,853 

0,716 

0,610 

0,526 

0,458 

0,403 

8500 

1,734 

1,370 

1,109 

0,917 

0,770 

0,657 

0,566 

0,493 

0,433 

9000 

1,866 

1,475 

1,195 

0,987 

0,830 

0,707 

0,609 

0,531 

0,467 

9500 

2,012 

1,589 

1,287 

1,064 

0,894 

0,762 

0,657 

0,572 

0,503 

10000 

2,170 

1,715 

1,389 

1,148 

0,965 

0,822 

0,709 

0,617 

0,543 

10500 

2,344 

1,852 

1,500 

1,240 

1,042 

0,888 

0,766 

0,667 

0,586 

11000 

2,323 

1,836 

1,487 

1,229 

1,033 

0,880 

0,759 

0,661 

0,581 

11500 

2,514 

1,986 

1,609 

1,330 

1,117 

0,952 

0,821 

0,715 

0,628 

12000 

2,720 

2,149 

1,741 

1,439 

1,209 

1,030 

0,888 

0,774 

0,680 

12500 

2,943 

2,325 

1,884 

1,557 

1,308 

1,115 

0,961 

0,837 

0,736 


Figure 5.6: Results lift coefficient 


The lift coefficient is influenced by density and true airspeed V (which is made up of Mach 
number and the speed of sound): 


C L 


2 ■ m ■ g 
p ■ M 2 ■ a 2 ■ S 


(5.5) 


This implies an increase in lift coefficient with an increase in altitude (because of the decreas¬ 
ing pressure). The second effect is that of true airspeed. A higher TAS results in a lower lift 
coefficient. A maximum TAS can be found at the highest Mach number and lowest altitude 
because of the increase of speed of sound with decreasing altitude. And in the same way a 
minimum TAS at the lowest Mach number and highest altitude. This effect can also be seen 
in Figure 5.6. An important thing to note is that there is no pursuit for maximum lift coeffi¬ 
cient during cruise flight, as this results in a high drag coefficient. The optimal lift coefficient 
values can be found in the yellow or orange zone on Figure 5.6. 
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Mach number 


0,4 

0,45 

0,5 

0,55 

0,6 

0,65 

0,7 

0,75 

0,8 


3000 

0,044 

0,033 

0,027 

0,023 

0,020 

0,019 

0,017 

0,017 

0,017 


3500 

0,048 

0,036 

0,029 

0,024 

0,021 

0,019 

0,018 

0,017 

0,018 


4000 

0,053 

0,039 

0,030 

0,025 

0,022 

0,020 

0,018 

0,018 

0,018 


4500 

0,058 

0,042 

0,033 

0,027 

0,023 

0,021 

0,019 

0,018 

0,018 


5000 

0,064 

0,046 

0,035 

0,029 

0,024 

0,022 

0,020 

0,019 

0,019 


5500 

0,071 

0,050 

0,038 

0,031 

0,026 

0,023 

0,021 

0,019 

0,019 


6000 

0,079 

0,055 

0,042 

0,033 

0,028 

0,024 

0,022 

0,020 

0,020 


6500 

0,088 

0,061 

0,046 

0,036 

0,030 

0,026 

0,023 

0,021 

0,021 


7000 

0,099 

0,068 

0,050 

0,039 

0,032 

0,027 

0,024 

0,022 

0,022 

0 ) 

“0 

7500 

0,112 

0,076 

0,056 

0,043 

0,035 

0,029 

0,026 

0,023 

0,023 

13 

■M 

8000 

0,127 

0,085 

0,062 

0,047 

0,038 

0,032 

0,028 

0,025 

0,024 

4-J 

< 

8500 

0,145 

0,097 

0,069 

0,052 

0,042 

0,035 

0,030 

0,027 

0,025 


9000 

0,165 

0,110 

0,078 

0,059 

0,046 

0,038 

0,032 

0,029 

0,027 


9500 

0,190 

0,125 

0,088 

0,066 

0,051 

0,042 

0,035 

0,031 

0,029 


10000 

0,218 

0,143 

0,100 

0,074 

0,057 

0,046 

0,039 

0,034 

0,031 


10500 

0,252 

0,164 

0,114 

0,084 

0,065 

0,052 

0,043 

0,037 

0,033 


11000 

0,248 

0,161 

0,112 

0,083 

0,064 

0,051 

0,043 

0,037 

0,033 


11500 

0,287 

0,186 

0,129 

0,094 

0,072 

0,058 

0,047 

0,040 

0,036 


12000 

0,334 

0,216 

0,149 

0,108 

0,082 

0,065 

0,053 

0,045 

0,040 

12500 

0,388 

0,250 

0,172 

0,124 

0,094 

0,074 

0,060 

0,050 

0,044 


Figure 5.7: Results drag coefficient 


The drag coefficient consists of three contributors: zero lift drag coefficient, wave drag in¬ 
crement, and induced drag coefficient. The induced drag coefficient is highly influenced by 
the lift coefficient as stated before. It is also the most important contributor, as can be seen on 
Figure 5.8. This figure shows the importance of each contributor that makes up the drag coef¬ 
ficient of an Airbus A320-200 at regular flight conditions (Mach 0.78 at 11 500 m altitude). 


The wave drag increment has a small influence on the total drag coefficient. It’s near zero at 
the lower Mach numbers and ramps up quickly when reaching the optimum cruise Mach 
number for this type of aircraft. But even when at maximum value, it makes up less than 10 % 
of the total drag coefficient. The only effect that is visible is that the drag coefficient stays at a 
constant value around the optimum cruise Mach number, instead of decreasing due to a de¬ 
creasing lift coefficient. 
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Induced drag; 
61,79% 



Wave drag; 5,33% 


Zero lift drag; 
32,88% 


Figure 5.8: Drag component contribution on the Airbus A320-200 

Due to the choice of method 1 for calculating the zero lift drag coefficient in the Excel tool 
(see Chapter 3.1.2), its calculation is mostly influenced by the Oswald factor. This means a 
decrease in C D0 when the Mach number increases, which corresponds with the theory; looking 
at the following equation (derived in Chapter 2), drag can be divided into two components 
which are influenced by speed in opposite ways; 



(2.3) 


These components are called parasitic drag and induced drag respectively. Parasitic drag de¬ 
pends on the zero lift drag coefficient and induced drag on the lift coefficient. Their behavior 
in function of speed is shown graphically in Figure 5.9: 
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Figure 5.9: Aircraft drag in function of flight speed (Spakovszky 2006) 


5.3.3 Fuel Consumption 

The following table and accompanying contour plot show the fuel consumption as kg of fuel 
per NM flown, per kg of aircraft MTOW. And this in function of Mach number and altitude. 


Mach number 



0,4 

0,45 

0,5 

0,55 

0,6 

0,65 

0,7 

0,75 

0,8 


3000 

8,05E-05 

7,26E-05 

6,95E-05 

6,94E-05 

7,14E-05 

7,59E-05 

8,23E-05 

9,16E-05 

l,09E-04 


3500 

8,28E-05 

7,35E-05 

6,95E-05 

6,86E-05 

6,99E-05 

7,36E-05 

7,92E-05 

8,76E-05 

l,03E-04 


4000 

8,54E-05 

7,48E-05 

6,98E-05 

6,80E-05 

6,86E-05 

7,16E-05 

7,65E-05 

8,40E-05 

9,84E-05 


4500 

8,85E-05 

7,64E-05 

7,03E-05 

6,78E-05 

6,77E-05 

7,00E-05 

7,41E-05 

8,07E-05 

9,38E-05 


5000 

9,21E-05 

7,85E-05 

7,13E-05 

6,79E-05 

6,70E-05 

6,86E-05 

7,20E-05 

7,79E-05 

8,97E-05 


5500 

9,62E-05 

8,09E-05 

7,25E-05 

6,82E-05 

6,67E-05 

6,76E-05 

7,03E-05 

7,53E-05 

8,60E-05 


6000 

l,01E-04 

8,38E-05 

7,42E-05 

6,90E-05 

6,66E-05 

6,68E-05 

6,89E-05 

7,32E-05 

8,28E-05 


6500 

l,06E-04 

8,71E-05 

7,62E-05 

7,00E-05 

6,69E-05 

6,64E-05 

6,78E-05 

7,13E-05 

7,99E-05 


7000 

l,12E-04 

9,10E-05 

7,87E-05 

7,15E-05 

6,76E-05 

6,63E-05 

6,70E-05 

6,99E-05 

7,74E-05 

<u 

"O 

7500 

l,19E-04 

9,54E-05 

8,16E-05 

7,33E-05 

6,86E-05 

6,66E-05 

6,66E-05 

6,88E-05 

7,54E-05 

=3 

■M 

8000 

l,26E-04 

l,00E-04 

8,50E-05 

7,56E-05 

6,99E-05 

6,72E-05 

6,66E-05 

6,80E-05 

7,37E-05 

'+-» 

< 

8500 

l,34E-04 

l,06E-04 

8,90E-05 

7,83E-05 

7,17E-05 

6,82E-05 

6,69E-05 

6,76E-05 

7,24E-05 


9000 

l,44E-04 

l,13E-04 

9,35E-05 

8,15E-05 

7,40E-05 

6,96E-05 

6,75E-05 

6,76E-05 

7,16E-05 


9500 

l,54E-04 

l,20E-04 

9,87E-05 

8,53E-05 

7,67E-05 

7,15E-05 

6,86E-05 

6,80E-05 

7,llE-05 


10000 

l,65E-04 

l,28E-04 

l,05E-04 

8,97E-05 

7,99E-05 

7,38E-05 

7,02E-05 

6,88E-05 

7,llE-05 


10500 

l,78E-04 

l,37E-04 

l,llE-04 

9,47E-05 

8,36E-05 

7,66E-05 

7,22E-05 

7,01E-05 

7,15E-05 


11000 

l,78E-04 

l,37E-04 

l,llE-04 

9,47E-05 

8,36E-05 

7,65E-05 

7,21E-05 

7,01E-05 

7,15E-05 


11500 

l,91E-04 

l,46E-04 

l,18E-04 

9,98E-05 

8,76E-05 

7,95E-05 

7,44E-05 

7,16E-05 

7,23E-05 


12000 

2,05E-04 

l,56E-04 

l,26E-04 

l,06E-04 

9,21E-05 

8,30E-05 

7,71E-05 

7,36E-05 

7,35E-05 


12500 

2,21E-04 

l,68E-04 

l,34E-04 

l,12E-04 

9,71E-05 

8,70E-05 

8,02E-05 

7,60E-05 

7,52E-05 


Figure 5.10: Results fuel consumption in kg fuel per NM flown, per kg of aircraft MTOW 
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Figure 5.11 : Contour plot of the fuel consumption results 


The color pattern of the aerodynamic efficiency results from Chapter 5.3.2 can clearly be rec¬ 
ognized in the fuel consumption results. The pattern is shifted towards higher Mach number 
because of the Mach number’s presence in the denominator of the fuel consumption equation 
(Eq. (5.3)). There is also a shift towards higher altitudes, caused by the decreasing TSFC 
(found in the nominator) with increasing altitude. 


c_g_ 

M-a-E 


fMTOW ~ 


(5.3) 
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Fuel consumption can also be represented in units of kg of fuel per 100 km flown, per seat. 



Mach number 


0,4 

0,45 

0,5 

0,55 

0,6 

0,65 

0,7 

0,75 

0,8 


3000 

2,13 

1,92 

1,84 

1,84 

1,89 

2,01 

2,18 

2,42 

2,88 


3500 

2,19 

1,94 

1,84 

1,81 

1,85 

1,95 

2,10 

2,32 

2,74 


4000 

2,26 

1,98 

1,85 

1,80 

1,82 

1,90 

2,02 

2,22 

2,60 


4500 

2,34 

2,02 

1,86 

1,79 

1,79 

1,85 

1,96 

2,14 

2,48 


5000 

2,44 

2,08 

1,89 

1,80 

1,77 

1,82 

1,91 

2,06 

2,37 


5500 

2,55 

2,14 

1,92 

1,81 

1,76 

1,79 

1,86 

1,99 

2,28 


6000 

2,67 

2,22 

1,96 

1,82 

1,76 

1,77 

1,82 

1,94 

2,19 


6500 

2,81 

2,30 

2,02 

1,85 

1,77 

1,76 

1,79 

1,89 

2,11 


7000 

2,96 

2,41 

2,08 

1,89 

1,79 

1,76 

1,77 

1,85 

2,05 

<u 

-Q 

7500 

3,14 

2,52 

2,16 

1,94 

1,81 

1,76 

1,76 

1,82 

1,99 

=3 

■M 

8000 

3,33 

2,66 

2,25 

2,00 

1,85 

1,78 

1,76 

1,80 

1,95 

■4—* 

< 

8500 

3,55 

2,81 

2,35 

2,07 

1,90 

1,80 

1,77 

1,79 

1,92 


9000 

3,80 

2,98 

2,47 

2,16 

1,96 

1,84 

1,79 

1,79 

1,89 


9500 

4,07 

3,17 

2,61 

2,26 

2,03 

1,89 

1,82 

1,80 

1,88 


10000 

4,38 

3,39 

2,77 

2,37 

2,11 

1,95 

1,86 

1,82 

1,88 


10500 

4,72 

3,63 

2,95 

2,50 

2,21 

2,03 

1,91 

1,85 

1,89 


11000 

4,72 

3,63 

2,94 

2,50 

2,21 

2,03 

1,91 

1,85 

1,89 

11500 

5,06 

3,87 

3,13 

2,64 

2,32 

2,10 

1,97 

1,89 

1,91 


12000 

5,43 

4,14 

3,33 

2,79 

2,44 

2,20 

2,04 

1,95 

1,94 


12500 

5,84 

4,43 

3,55 

2,96 

2,57 

2,30 

2,12 

2,01 

1,99 


Figure 5.12: Results fuel consumption in kg fuel per 100 km flown, per seat 


As a comparison to these results, an average car can be said to consume about 6 1 of gasoline 
per 100 km driven on the highway. With an average density of gasoline of about 0.755 kg/1, 
this equals 4.53 kg of gasoline. This means that one person driving in a car would consume 
more fuel than one person traveling on this aircraft at optimal flight conditions or even regular 
fight conditions (Mach 0.78 at 11 500 m). Keep in mind that this is greatly simplified: the 
LTO cycle is not accounted, an aircraft is used for longer ranges than a car in general, 
transport to and from the airport is not accounted, etc. 


5.4 Conclusion Fuel Consumption 

In this chapter fuel consumption is calculated by use of aerodynamic parameters and the air¬ 
craft’s TSFC. With the aerodynamic parameters based on flight parameters and (a minimum 
of) geometric parameters of the aircraft. The TSFC is based on flight parameters and basic 
engine parameters. All results are for an Airbus A320-200, but the methods used are meant 
for any passenger jet aircraft. 

All models used for calculating aerodynamic parameters have one or more factors based on 
statistical data. This introduces a certain level of inaccuracy in the values of individual results 
(e.g. one result for a certain Mach number of altitude). However, the scope of this work is not 
to generate exact data for fuel consumption, but to visualize the trend of fuel consumption 
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variation when altering flight parameters. This goal is still accomplished because the same 
models and statistical factors are used for each individual result, ensuring that the fuel con¬ 
sumption variation stays credible. In addition to this, TSFC results seem to be within bounds 
of actual engine data (see Chapter 5.3.1) and the fuel consumption results are within the ex¬ 
pected order of magnitude for a single aisle aircraft. 

The fuel consumption results in Chapter 5.3.3 show a big influence of the aerodynamic effi¬ 
ciency E. They follow the colored shape of the optimal ratio of lift coefficient to drag coeffi¬ 
cient. Aerodynamic efficiency can be optimized by increasing the lift coefficient while reduc¬ 
ing the drag coefficient. The TSFC, Mach number, and local speed of sound have a smaller 
influence on fuel consumption. 

The optimal flight parameters for minimal fuel consumption of an Airbus A320-200, calculat¬ 
ed with the methods present, are Mach 0.65 at an altitude of 7000 m. However, the whole 
dark green zone in Figure 5.10 only differs from this minimum with up to 2.45 %. When 
looking at the regular flight parameters of Mach 0.78 at 11 500 m, the fuel consumption dif¬ 
fers 7.01 % from the minimum. 

The importance of considering fuel consumption in function of Mach number and altitude is 
that both flight parameters can be changed without changing anything to the aircraft’s design. 
At least for small variations in cruise speed or altitude. When flying at greatly reduced alti¬ 
tudes in comparison to common jet cruise altitudes, the wing loading of the aircraft could be 
changed to further reduce fuel consumption. This scenario is discussed in Chapter 8.3. 
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6 Flight Parameter Optimization: DOC 

6.1 General 

The direct operating costs represent the economic impact of operating an aircraft. In this chap¬ 
ter DOC will be represented as an annual cost and as the seat-mile cost. DOC is composed of 
seven categories and several subcategories as explained in Chapter 3.3. The aim is to visualize 
the economic impact of flying at varying flight parameters. To have a generic model and 
maintain the option to compare different aircraft, the seat-mile cost will be used as the prima¬ 
ry metric. 


6.2 Important Parameters 

The most influential parameters in DOC calculation are represented in Figure 6.1. Inputs are 
green, important intermediate parameters in blue and the final results in grey. DOC is divided 
into time-independent and time-dependent costs (see Chapter 3.3). The time-independent 
costs in blue on the left are only influenced by the type of aircraft, and airframe- and engine 
prices. The type of aircraft detennines the depreciation period and interest rate. The time- 
dependent costs are influenced by flight time, fuel price, crew and maintenance costs, and dif¬ 
ferent fees and charges. Only flight time and fuel price are actually varying. Crew and 
maintenance costs and fees are fixed 1989 values which are adjusted for inflation. 
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Figure 6.1 : Flowchart of DOC calculation 


6.2.1 DOC Range 

The DOC range R DO c is considered as a typical range for the type of aircraft considered. It is 
defined as 50 % of the range for maximum payload in Scholz 2015. The range for maximum 
payload can be found on the aircraft’s payload-range diagram, which is publicly available in 
the aircraft’s “Documents for airport planning” (see Airbus 2005a for the A320 document). 
Range has a big influence on DOC in two main ways: it influences flight time and it influ¬ 
ences fuel consumption. Its influence on flight time is discussed in Chapter 6.2.2. The influ¬ 
ence on fuel consumption can be visualized with the so-called bath tub curve (Figure 6.2). 

This curve displays the relationship between fuel consumption and range. At short range, fuel 
consumption is relatively high due to the necessary amount of reserve fuel and alternate fuel 
that must be on board independent of range. With increasing range, fuel consumption de¬ 
creases up to an optimum range and starts increasing again slightly after this point. From the 
optimum point on, the aircraft must take on more fuel in order to transport the fuel required 
for the longer range. At a certain point fuel consumption increases in a steep way due to pay- 
load restrictions. At this point the aircraft must reduce payload (passengers or cargo) in order 
to stay under the maximum take-off weight (Burzlaff 2017). 
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Figure 6.2: Bath tub curve of an exemplary aircraft (Burzlaff 2017) 


6.2.2 Flight Time 

Flight time tf is the most important parameter when considering seat-mile cost. The flight 
time in hours could be considered as only dependent on true airspeed and DOC range and 
could be calculated as follows: 


Rdoc 


( 6 . 1 ) 


However, this equation is a simplification. It calculates the flight time as if the whole range 
would be flown at cruise speed. A more correct approach is to include the speed differences 
during the climb and descend flight stages (take-off and landing only take up a very small 
amount of time compared to the total flight time and will be neglected). Especially for short 
ranges climb and descent time takes up a considerable share of the total flight time. The added 
flight time can be calculated from the climb rate and descend rate at certain altitude zones 
with the corresponding speeds. These rates and speeds are aircraft-dependent and regulated by 
the authorities. 


Two ways to calculate flight time are discussed here. The first one is to look for aircraft spe¬ 
cific data in the aircraft performance database by Eurocontrol 2019 or their (less up-to-date) 
BADA tables (Eurocontrol 1998) and calculate the flight time per colored segment. See Fig¬ 
ure 6.3 for the Airbus A320-200 performance data. 
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MACH CLIMB 


ROC: lOOOft/min 
CLIMB (to FL240) 

IAS: 290 kts 
ROC: 1400 ft/min 
CLIMB (to FL150) 

IAS: 290 kts 
ROC 2000 ft/min 



INITIAL CLIMB 
(to 5000 ft) 


IAS: 175 kts 
ROC: 2500 ft/min 


TAKE-OFF 
MTOW 73900 kg 

WTC: M / RECAT-EU: Upper Medium 

Distance: 21 90 m 

V2 (IAS): 145 kts 




to FL24 



LANDING 


Vat(IAS): 137 kt 
A PC: c 
Distances 440 m 


Figure 6.3: Aircraft performance for the Airbus A320-200 (Eurocontrol 2019) 


The second option is to use a more generic approach and work with common practice and 
aviation regulations as a guideline to fonn an averaged model for approximating flight time of 
any given passenger jet aircraft. Assumptions regarding aircraft perfonnance data are based 
on a comparison of BAD A tables of four popular aircraft (Eurocontrol 1998): Airbus A320- 
200, Airbus A340, Boeing 737-200, and Boeing 747-400. Corresponding graphs can be found 
in Appendix D. 


The following assumptions are made when estimating flight time (full calculation can be 
found in the ‘Flight time’ sheet in the Excel tool): 


- LTO cycle is neglected 

- Maximum IAS of 250 knots below FL100 by regulation 

- Average ROC and ROD of 2 000 ft/min 

- TAS varies linearly between 150 knots and 290 knots under FL100 

- Between FF120 and FF300: variation between 360 knots and cruise speed 

- Above FF300: constant cruise speed 


The difference in DOC results by using this model instead of Eq. (6.1) is given in Figure 6.4. 
All calculations are done for the Airbus A320-200 at an altitude of 11 000 m. A difference of 
up to 7.9 % can be seen in the actual flight time at higher Mach numbers. This is because of 
the increased relative duration of the climbing and descending stages at shorter flight times. 
The annual DOC varies less than 3.5 % between the two methods, this could be considered 
acceptable. The seat-mile cost, which is used in the final combined model (see Chapter 8), 
varies up to about 4.1 %. This could be considered within boundaries given the not so accu- 
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rate nature of the DOC model. But in an effort to increase the correctness of the DOC or seat- 
mile cost model, the flight time estimation model will be used instead of Eq. (6.1). 



Figure 6.4: Difference in percentage between flight time estimation model and simple equation 


6.2.3 Number of Flights per Year 

The number of flights per year is an important parameter in calculating the annual DOC cost. 
It is used in the calculation of fuel cost, maintenance cost, staff cost, and fees and charges. Es¬ 
sentially every cost that is first calculated per flight and later changed to an annual cost. It is 
also a necessary factor to reform the annual DOC cost into seat-mile cost. The number of 
flights per year is based on annual aircraft utilization and flight time. And in turn the annual 
aircraft utilization U a j is also based on flight time: 


L u 1 


Ua ’ f tf 't f + k 


U 2 


The number of flights per year n t a can then be defined as: 


U a ,f _ 

£/ tf + k U2 


( 6 . 2 ) 


(6.3) 


With k ai and k U2 based on the type of aircraft: 
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Table 6.1: Factors for aircraft utilization equation (Scholz 2015) 


Aircraft type 

kUl (in h) 

kU2 (in h) 

Short-/medium-range 

3750 

0,75 

Long-range 

4800 

0,42, 


Annual aircraft utilization and number of flights per year can be reshaped into daily values by 
dividing both by 365. Daily aircraft utilization and number of flights per day can be used to 
compare aircraft or airlines and, in this case, to verify the correctness of the user’s results 
when calculating annual aircraft utilization in the Excel tool. Data for American Airlines from 
Swelbar 2018 indicates the following averages: 


Table 6.2: Average daily aircraft utilization (Swelbar 2018) 


Aircraft type 

Daily aircraft utilization (h) 

Number of flights per day 

Short-/medium-range 

8,2 

4,5 

Long-range 

11,5 

1,5 


6.2.4 Seat-Mile Cost 


The seat-mile cost is the direct operating cost for one seat in the aircraft over a distance of one 
NM. This metric is often used in an economic context to compare different aircraft models, 
independent of their fuselage size or design range. It is calculated as follows: 


r — 

u s,m 


c, 


a/c,a 


n 


seat 


n 


t,a 


R 


DOC 


(3.11) 


6.2.5 Dollar Prices and Inflation 

Every currency value in the DOC model is represented in U.S. Dollar (USD). All estimated 
USD values in the AEA1989 model are based on 1989 values. This means a correction is 
needed to achieve present date USD values. A compensation for inflation must be used. The 
inflation factor is defined as: 


k INF = (1 + ViN F ) nyear nmethod (6.4) 

The annual mean inflation rate p INF for the period 1989 - 2017 in the United States is 
2.5547 % (Worldbank 2019). Multiplying all USD values in the model with an inflation fac¬ 
tor is a basic rough compensation. In reality prices change due to a multitude of effects. Un¬ 
ions can demand higher crew wages; airport fees can change, or aircraft manufacturers can al¬ 
ter their prices. This means that the compensation for inflation only gives an approximate val- 













68 


ue. Fuel prices for example can hardly be approximated by an inflation factor, due to oil crises 
or political events strongly influencing the oil price. For this reason, the current fuel price 
must be added by the user. In conclusion, this correction for inflation means that all USD val¬ 
ues in the DOC model are only an approximation. 


6.3 Results DOC 

DOC results can be represented as annual DOC or as seat-mile cost. Annual DOC returns a 
clear number of the cost paid by the airline for operating the aircraft, but this value differs a 
lot between aircraft of different sizes. Seat-mile cost is therefore the preferred output to keep 
the model as generic as possible. By dividing annual DOC by number of flights per year, 
DOC range and the number of seats for the seat-mile cost, different types of aircraft can be 
compared independent of their size. 


6.3.1 Flight Time 


The flight time is shown in hours. 



Mach number 


0,4 

0,45 

0,5 

0,55 

0,6 

0,65 

0,7 

0,75 

0,8 


3000 

3,17 

2,82 

2,55 

2,33 

2,15 

2,00 

1,87 

1,75 

1,65 


3500 

3,22 

2,87 

2,60 

2,38 

2,20 

2,04 

1,91 

1,80 

1,70 


4000 

3,24 

2,89 

2,62 

2,40 

2,22 

2,06 

1,93 

1,82 

1,72 

4500 

3,26 

2,90 

2,63 

2,41 

2,23 

2,08 

1,95 

1,83 

1,73 


5000 

3,28 

2,92 

2,65 

2,43 

2,25 

2,09 

1,96 

1,85 

1,75 


5500 

3,30 

2,94 

2,66 

2,45 

2,26 

2,11 

1,98 

1,86 

1,76 


6000 

3,33 

2,96 

2,68 

2,46 

2,28 

2,12 

1,99 

1,88 

1,77 

6500 

3,35 

2,98 

2,70 

2,48 

2,29 

2,14 

2,00 

1,89 

1,79 

E 

7000 

3,37 

3,00 

2,71 

2,49 

2,30 

2,15 

2,01 

1,90 

1,79 

(U 

7500 

3,39 

3,02 

2,73 

2,50 

2,32 

2,16 

2,02 

1,90 

1,80 

=5 

8000 

3,42 

3,04 

2,75 

2,52 

2,33 

2,17 

2,03 

1,91 

1,81 

< 

8500 

3,44 

3,06 

2,76 

2,53 

2,34 

2,18 

2,04 

1,92 

1,81 


9000 

3,46 

3,08 

2,78 

2,55 

2,35 

2,19 

2,05 

1,92 

1,82 


9500 

3,49 

3,10 

2,80 

2,58 

2,39 

2,23 

2,10 

1,98 

1,88 


10000 

3,51 

3,13 

2,82 

2,59 

2,41 

2,25 

2,11 

1,99 

1,89 


10500 

3,54 

3,15 

2,84 

2,61 

2,42 

2,26 

2,13 

2,01 

1,90 


11000 

3,56 

3,17 

2,86 

2,63 

2,44 

2,28 

2,14 

2,02 

1,91 


11500 

3,56 

3,17 

2,86 

2,63 

2,44 

2,28 

2,14 

2,02 

1,92 


12000 

3,56 

3,17 

2,86 

2,63 

2,44 

2,28 

2,14 

2,02 

1,91 


12500 

3,56 

3,17 

2,86 

2,63 

2,44 

2,28 

2,14 

2,02 

1,92 


Figure 6.5: Results flight time in hours 


The flight time is a result of TAS and cruise altitude. The dependency on cruise altitude (or 
flight level) can be seen in Figure 6.6. This means that different altitude zones will result in 
bigger or smaller variations in flight time, as can be seen for example when varying cruise al- 
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titude from 3000 m to 3500 m (a shift from around FL100 to FL115). This results in a sub¬ 
stantial difference in flight time compared to other altitude zones. A lower cruise altitude also 
implicates less time spent during climb and descent. 

The TAS itself depends on Mach number and altitude. A higher Mach number logically re¬ 
flects a higher TAS. And the altitude influences the local speed of sound, which in turn influ¬ 
ences the TAS. A higher altitude means a lower TAS for the same Mach number. 



Figure 6.6: Airbus A320-200 TAS during climb and descent (approximation) 
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6.3.2 Annual DOC 


The annual DOC is represented as a yearly cost in USD for the operation of one aircraft. 



Mach number 


0,4 

0,45 

0,5 

0,55 

0,6 

0,65 

0,7 

0,75 

0,8 


3000 

3,10E+07 

3,17E+07 

3,24E+07 

3,31E+07 

3,39E+07 

3,47E+07 

3,55E+07 

3,64E+07 

3,75E+07 


3500 

3,09E+07 

3,16E+07 

3,22E+07 

3,29E+07 

3,36E+07 

3,44E+07 

3,52E+07 

3,60E+07 

3,71E+07 


4000 

3,09E+07 

3,15E+07 

3,22E+07 

3,29E+07 

3,35E+07 

3,43 E+07 

3,50E+07 

3,58E+07 

3,68E+07 


4500 

3,09E+07 

3,15E+07 

3,22E+07 

3,28E+07 

3,35E+07 

3,41E+07 

3,48E+07 

3,56E+07 

3,66E+07 


5000 

3,09E+07 

3,15E+07 

3,21E+07 

3,27E+07 

3,34E+07 

3,40E+07 

3,47E+07 

3,54E+07 

3,63E+07 


5500 

3,09E+07 

3,15E+07 

3,21E+07 

3,27E+07 

3,33E+07 

3,39E+07 

3,46E+07 

3,53E+07 

3,62E+07 


6000 

3,10E+07 

3,15E+07 

3,21E+07 

3,27E+07 

3,33E+07 

3,39E+07 

3,45E+07 

3,52E+07 

3,60E+07 


6500 

3,10E+07 

3,15E+07 

3,21E+07 

3,26E+07 

3,32E+07 

3,38E+07 

3,44E+07 

3,51E+07 

3,58E+07 

E 

7000 

3,llE+07 

3,16E+07 

3,21E+07 

3,26E+07 

3,32E+07 

3,38E+07 

3,43E+07 

3,50E+07 

3,57E+07 

(U 

7500 

3,12E+07 

3,16E+07 

3,21E+07 

3,26E+07 

3,32E+07 

3,37E+07 

3,43E+07 

3,49E+07 

3,56E+07 

13 

■M 

8000 

3,12E+07 

3,17E+07 

3,21E+07 

3,26E+07 

3,31E+07 

3,37E+07 

3,43E+07 

3,48E+07 

3,55E+07 

< 

8500 

3,12E+07 

3,17E+07 

3,22E+07 

3,26E+07 

3,31E+07 

3,37E+07 

3,42E+07 

3,48E+07 

3,55E+07 


9000 

3,12E+07 

3,18E+07 

3,22E+07 

3,27E+07 

3,32E+07 

3,37E+07 

3,42E+07 

3,48E+07 

3,54E+07 


9500 

3,llE+07 

3,19E+07 

3,23E+07 

3,27E+07 

3,31E+07 

3,35E+07 

3,40E+07 

3,45E+07 

3,51E+07 


10000 

3,llE+07 

3,19E+07 

3,24E+07 

3,27E+07 

3,31E+07 

3,35E+07 

3,40E+07 

3,45E+07 

3,50E+07 


10500 

3,llE+07 

3,19E+07 

3,24E+07 

3,28E+07 

3,32E+07 

3,36E+07 

3,40E+07 

3,44E+07 

3,50E+07 


11000 

3,10E+07 

3,18E+07 

3,24E+07 

3,27E+07 

3,31E+07 

3,35E+07 

3,39E+07 

3,44E+07 

3,49E+07 


11500 

3,llE+07 

3,18E+07 

3,25E+07 

3,28E+07 

3,32E+07 

3,36E+07 

3,40E+07 

3,44E+07 

3,49E+07 


12000 

3,llE+07 

3,19E+07 

3,26E+07 

3,30E+07 

3,33E+07 

3,37E+07 

3,41E+07 

3,45E+07 

3,50E+07 


12500 

3,llE+07 

3,19E+07 

3,26E+07 

3,31E+07 

3,34E+07 

3,38E+07 

3,42E+07 

3,46E+07 

3,50E+07 


Figure 6.7: Results annual DOC for one aircraft in USD 


The annual DOC results follow the inverse pattern of the flight time, because a higher flight 
time means a lower number of flights per year. The number of flights per year has a heavy in¬ 
fluence on all time-dependent (or cost index) components of the DOC: fuel, maintenance, 
crew costs, and fees and charges. All of these costs decrease with a lower number of flights 
per year. 


Another observation is that the annual DOC follows the pattern of minimal fuel consumption 
slightly. Costs reach their lowest value in each column around the point of minimum fuel con¬ 
sumption for that column. The strength of this relationship depends on the actual fuel price 
and flight mission of the aircraft. A higher fuel price shows a clearer pattern towards mini¬ 
mum fuel consumption in the annual DOC results. Furthermore, short-/medium-range aircraft 
like the Airbus A320-200 depend less on fuel consumption than long-range aircraft. In the re¬ 
sult for the Airbus A380 in Appendix E, a bigger variation is visible in annual DOC results 
per single column. 


The annual DOC result conveys the idea that flying lower and slower would be an optimal 
flight condition. There are however two reasons that this is not the case. The first one is that 
flying slower would increase the flight time. At Mach 0.4 this would be almost double the 
flight time of a passenger jet aircraft nowadays. Passenger willingness to buy tickets for such 
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a flight would be substantially lower (this idea is further elaborated in Chapter 6.4). A second 
reason is that the increased flight time leads to a lower number of flights per year. This in turn 
negatively impacts the seat-mile cost and the ability of the airline to make a profit by selling 
tickets. This is one of the reasons why seat-mile cots is chosen as the representative economic 
metric instead of annual DOC. 

The distribution of all DOC components is shown in Figure 6.8 for the regular flight condi¬ 
tions of the Airbus A320-200 (cruise speed of 0.78 M at 11 500 m altitude). 



Maintenance 

15% 

Figure 6.8: Annual DOC distribution for Mach 0.78 at 11 500 m 

The first thing that stands out in the pie chart is the big slice of fees and charges. This seems 
abnormal but was researched by Meyer 2004 and deemed to be realistic. The fees and charges 
share consists of landing fees, navigation or ATC fees, and ground handling fees. All USD 
values here are corrected for inflation from the 1989 values provided in the AEA1989 meth¬ 
od. 

A second observation is the low contribution of fuel costs to the total DOC. This is a result of 
the current low price of jet fuel (around 1.8 USD per gallon at the time of publication) and the 
short DOC range for this aircraft. See Figure 2.4 for the jet fuel price evolution over the last 
decades. Another reason for the low fuel cost is the fact that this is a short-/medium-range air¬ 
craft. The fuel cost share increases significantly once the DOC range is increased (see Appen¬ 
dix E for Airbus A380 case). 

Figure 6.9 presents an overview of annual DOC distribution for four different flight parameter 
scenarios. The first observation is that depreciation, interest, and insurance stay at a constant 
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value since they only depend on the type of aircraft (short-/medium-range or long-range) and 
not on the flight time. Furthermore, fees and charges increase substantially with increasing 
Mach number due to the increase in number of flights per year (shorter flight time). This im¬ 
plies a greater number of LTO cycles and therefore an increase in all fees and charges: land¬ 
ing, navigation, and ground handling. 


Fuel costs increase with number of flights per year and increasing fuel consumption to a lesser 
extent. Maintenance costs follow the same increase with an increasing number of flights per 
year. Staff costs, in contrary, decrease slightly due to the lower number of block hours. 



Case 1 Case 2 

Case 3 


Case 4 

M 

0,4 

0,5 

0,75 

0,8 

H (no) 

3000 

7000 

11000 

12500, 

Depreciation 

$6 200 021,45 

$6 200 021,45 

$6 200 021,45 

$6 200 021,45 

Interest 

$5 101 928,76 

$5 101 928,76 

$5 101 928,76 

$5 101 928,76 

Insurance 

$423 617,79 

$423 617,79 

$423 617,79 

$423 617,79 

Fuel 

$2 413 909,79 

$2 691 784,30 

$3 128 502,36 

$3 396 793,16 

Maintenance 

$5 023 752,69 

$5 106 346,31 

$5 285 011,84 

$5 320 396,10 

Staff 

$5 130 719,38 

$5 031 588,64 

$4 819 748,91 

$4 778 910,89 

Fees and charges 

$6 674 835,27 

$7 553 653,71 

$9 443 810,82 

$9 813 472,54 


Flight time (h) 

3,17 

2,71 

2,02 

1,92 

Seat-mile costs 

$0,266 

$0,244 

$0,209 

$0,205 


$ 40,00 


$ 35,00 


$ 30,00 


$ 25,00 


$ 20,00 


$ 15,00 


$ 10,00 


$ 5,00 


$ 0,00 





__ 







$ 6,67 


$ 7,55 


$ 9,44 





$ 5,13 


$ 5,03 


$ 4,82 


$ 4,78 



$ 5,02 


$ 5,11 


$ 5,29 


$ 5,32 



$ 2,41 


$ 2,69 


$ 3,13 

Crwin 


$ 3,40 







$ 5,10 


$ 5,10 


$ 5,10 


$ 5,10 



$ 6,20 


$ 6,20 


$ 6,20 


$ 6,20 



12 3 

■ Depreciation ■ Interest ■ Insurance ■ Fuel ■ Maintenance ■ Staff ■ Fees and charges 


Figure 6.9: 


Breakdown of annual DOC for four flight parameter scenarios 
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6.3.3 Seat-Mile Cost 

The seat-mile cost is shown as the cost in USD for one seat in the aircraft flown over a dis¬ 
tance of one NM. 


Mach number 



0,4 

0,45 

0,5 

0,55 

0,6 

0,65 

0,7 

0,75 

0,8 


3000 

0,266 

0,248 

0,235 

0,224 

0,216 

0,209 

0,204 

0,200 

0,198 


3500 

0,269 

0,251 

0,237 

0,226 

0,218 

0,211 

0,205 

0,201 

0,199 


4000 

0,271 

0,252 

0,238 

0,227 

0,218 

0,211 

0,206 

0,202 

0,199 


4500 

0,272 

0,253 

0,239 

0,228 

0,219 

0,212 

0,206 

0,202 

0,199 


5000 

0,274 

0,254 

0,240 

0,229 

0,220 

0,213 

0,207 

0,202 

0,199 


5500 

0,276 

0,256 

0,241 

0,229 

0,220 

0,213 

0,207 

0,203 

0,199 


6000 

0,277 

0,257 

0,242 

0,230 

0,221 

0,214 

0,208 

0,203 

0,199 


6500 

0,279 

0,258 

0,243 

0,231 

0,222 

0,214 

0,208 

0,203 

0,200 

E 

7000 

0,281 

0,260 

0,244 

0,232 

0,222 

0,215 

0,208 

0,203 

0,199 

(U 

7500 

0,283 

0,262 

0,245 

0,233 

0,223 

0,215 

0,209 

0,203 

0,200 

=5 

+-> 

8000 

0,286 

0,264 

0,247 

0,234 

0,224 

0,216 

0,209 

0,204 

0,200 

< 

8500 

0,287 

0,265 

0,248 

0,235 

0,225 

0,216 

0,209 

0,204 

0,200 


9000 

0,288 

0,268 

0,250 

0,237 

0,226 

0,217 

0,210 

0,204 

0,200 


9500 

0,290 

0,270 

0,252 

0,239 

0,228 

0,220 

0,213 

0,207 

0,203 


10000 

0,291 

0,272 

0,253 

0,240 

0,229 

0,221 

0,214 

0,208 

0,203 


10500 

0,293 

0,273 

0,255 

0,242 

0,231 

0,222 

0,215 

0,209 

0,204 


11000 

0,294 

0,274 

0,256 

0,243 

0,232 

0,223 

0,215 

0,209 

0,204 


11500 

0,294 

0,274 

0,258 

0,243 

0,232 

0,223 

0,216 

0,209 

0,205 


12000 

0,294 

0,274 

0,258 

0,244 

0,233 

0,224 

0,216 

0,210 

0,205 


12500 

0,295 

0,274 

0,259 

0,245 

0,234 

0,224 

0,217 

0,210 

0,205 


Figure 6.10: Results seat-mile cost in USD 





Figure 6.11 : Contour plot of the seat-mile cost results 
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For seat-mile cost a general trend is visible towards high Mach number and low altitude. This 
follows the pattern of flight time entirely. Lower altitude means higher TAS for the same 
Mach number. Flight time is present in the seat-mile cost equation in the fonn of number of 
flights per year. Annual DOC evolves negatively towards lower altitude, but the influence of 
the number of flights per year on the seat-mile cost is higher and thus nullifies the negative 
influence of annual DOC. 
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When looking at the seat-mile cost equation, annual DOC and number of trips per year are the 
only variable parameters when considering a single aircraft. DOC range could be altered to 
influence the seat-mile cost, but this is a question of operational flight planning done by the 
airline. 


6.4 Added Value 

In the annual DOC results can be seen that flying low and slow results in the lowest annual 
direct operating costs. DOC is purely built up out of economic data and questions can be 
asked here about the human factor in these results. When considering that a passenger jet air¬ 
craft would fly at about half the speed normally flown, would passengers still be willing to 
pay the same ticket price for a flight that takes twice as long? This customer willingness is 
difficult to put into numbers because it is highly subjective. In the work of Nita 2013 a tenn 
called added value is used. This is an effort to add non-economic factors like performance, 
operating flexibility, or passenger comfort to the economic equation. 

To decide on the weight of each of the added value factors, 22 questionnaires were taken by 
Nita from experts, aircraft design PhD students and aircraft design students. The final weights 
were determined from the best answer (i.e. best consistency and high coefficients of determi¬ 
nation) corrected by technical insight, expert views and the average of all answers (Nita 
2013). An overview of the added value factors and their final weights is given in Table 6.3. 
Cruise speed is the largest contributing factor to the added value, followed by passenger com¬ 
fort factors like seat pitch and -width or overhead bin volume. 
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Table 6.3: Final attributed weights to the added value factor (Nita 2013) 


Economics 


Added 

Values 


75 


25 


Absolute Added 
weights Values 


% 


Equiv. ton- 
mile costs 

100 




% 




% 

scaled to 

75.00% 100 % 



Performance 

35 



Airport performance 

50 

% 


Landing field length 

0 

% 

0.00% 0.00% 

3.50% 14.00% 

0.88% 3.50% 

4.38% 17.50% 

3.30% 13.20% 

2.20% 8.80% 
1.10% 4.40% 

0.55% 2.20% 

0.55% 2.20% 

2.20% 8.80% 
1.10% 4.40% 

Take-off field length 

80 

% 

Relative landing weight (nWiriMTo) 

20 

% 

Cruise performance 

50 

% 


Cruise speed 100 

% 

Passenger 

Comfort 

55 



Concerning all 
passengers 

80 



Seat pitch 

30 

% 

Seat width 

20 

% 

Armrest width 

10 

% 

Aisle width 

5 

5 

% 

Aisle height 

% 

Overhead bin volume per pax 

20 

% 

Aircraft gust sensibility 

10 

% 

Concerning part of 

20 


Sidewall clearance 

10 % 

0.28% 1.10% 







the passengers 



Number of "excuse-me" seats 

90 % 

2.48% 9.90% 

2.00% 8.00% 

0.25% 1.00% 

0.25% 1.00% 

Cargo 

Handling 

10 

% 


Concerning cargo 

80 



Containerized cargo (yes/no) 

too 

% 

Concerning cargo 
working conditions 

20 



Accessibility factor 

Cargo compartment height 

s 

1 - 


Check: 100.00% 100.00% 


Although the added value approach does not give a comprehensive answer to the question if 
passengers would still be willing to pay the same ticket price for a flight that takes twice as 
long, it shows that the cruise speed of an aircraft is considered as an important added value. In 
conclusion, this is another reason why seat-mile cost is a more appropriate metric for calculat¬ 
ing optimum economic benefit. 


6.5 Conclusion DOC 

In this chapter the economic impact of flying at different flight parameters is analyzed. The 
full method is based on the DOC method of Scholz 2015, which is in turn based on the 
AEA1989 method. The results are expressed in the form of annual cost and seat-mile cost of 
one aircraft. 

The accuracy of the model largely depends on 1989 USD values for prices, fees, and charges 
adjusted for an average inflation or depends on statistical data. This produces a certain inaccu¬ 
racy in the actual USD values. But again, as in Chapter 5, this does not pose a problem when 
only comparing different USD values for varying flight parameters to find a trend in DOC 
variation. 

In general, the annual DOC and seat-mile cost follow a different and almost inverse pattern. 
Annual DOC decreases when flying slower and has a minimum value for each Mach number 
at around the point of minimum fuel consumption. The minimum value would be even more 
emphasized by higher fuel prices or for a long-range aircraft. This result conveys the idea that 
flying at a lower Mach number is beneficial in an economic context. However, this is not the 
case due to several reasons: a low passenger-willingness to pay for very slow flights (see 
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Chapter 6.4 about added value), a cutback in number of flights per year and the corresponding 
ticket income for the airline, and finally an increase in seat-mile cost for the same aircraft. 

The main reason the annual DOC decreases with longer flight times is due to the lower num¬ 
ber of flights per year that are possible and the corresponding lower fuel costs, maintenance 
costs, staff costs, and fees and charges. This is a logical consequence: operating any vehicle to 
transport passengers less than its full potential (be it an aircraft or some other means of 
transport) will decrease its operational costs. But this defeats the main purpose of the vehicle: 
generating as much revenue as possible during its limited lifespan. 

Seat-mile cost follows the pattern of flight time. A minimal flight time results in a minimal 
seat-mile cost. This is because a lower flight time enables a higher number of flights per year. 
And therefore, cruise speed or cruise Mach number has a big influence on seat-mile cost. At 
an altitude of 11 500 m, the seat-mile cost is 30.47 % lower when flying at Mach 0.8 com¬ 
pared to Mach 0.4. Compared to this, the influence of cruise altitude is less distinct. It varies 
between 3.41 % and 9.58 % from high to low Mach numbers respectively. 

It is clear that seat-mile cost portrays the economic effect of varying flight parameters better 
than the annual DOC. Imagine a fixed ticket price for a certain flight mission. A lower seat- 
mile cost means a higher return on this ticket. In this way, the lowest flight time ensures the 
optimal economic benefit for an aircraft. This results in a general outcome for DOC: a higher 
true airspeed during cruise (higher Mach number and lower altitude) results in the greatest 
number of flights per year and consequently in the lowest seat-mile cost possible. An added 
benefit of high TAS being at the lowest cruise altitude is that climb- and descent stages take 
up less time in this case. Which in turn further improves flight time. 

Another observation that can be made from the results in this chapter is that fuel consumption 
has a minimal influence on annual DOC and especially on the seat-mile cost for this type of 
aircraft. As mentioned before, this is as a result of the current relatively low fuel price and be¬ 
cause the Airbus A320-200 is a short-/medium-range aircraft. Meaning that it typically spends 
more time on the ground loading and unloading and less time in the air than a typical wide- 
body long-range aircraft (Swelbar 2018). Airlines could therefore be inclined towards less ef¬ 
ficient (and thus probably less expensive) aircraft over expensive and efficient aircraft for 
their short- to medium-range flights. 
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7 Flight Parameter Optimization: Environmental 
Impact 

7.1 General 

The environmental impact of an aircraft can be defined in a number of different ways. The 
model is already narrowed down by the fact that only cruise flight is considered. Local air 
quality and noise pollution are not applicable at cruise altitudes. For this model, environmen¬ 
tal impact is considered as a combination of resource depletion and engine emissions. With 
resource depletion equal to fuel consumption. Engine emissions are estimated with equivalent 
CO 2 mass, a metric used to combine the effects of CO 2 , NO x and AIC into one equation. 
Equivalent CO 2 mass is introduced and discussed in Chapter 3.5. 


7.2 Important Parameters 

The most influential parameters in the environmental impact calculation are represented in 
Figure 7.1. Inputs are green, important intermediate parameters in blue and the final results in 
grey. CO 2 and NO x emissions are influenced by the aircraft’s fuel consumption, AIC emission 
only by its range. 



Figure 7.1: 


Flowchart of environmental impact calculation 
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7.2.1 Emission Index 

The emission index of a specie is the mass (in g or kg) of this specie emitted per kg of fuel 
burned. The emission indices that are considered here are EICO 2 and EI_NO x . The emission 
index of aircraft-induced clouds is hard to define because the formation of contrails depends 
on the altitude and distance flown and not on the mass of fuel burned. The emission index for 
CO 2 is independent of cruise altitude and equals a constant 3.16 kg CO 2 per kg of fuel. The 
emission index for NO x on the other hand is dependent on the cruise altitude and is calculated 
with the Boeing fuel flow method 2 (Baughcum 1996). Application of this method in the Ex¬ 
cel tool can be found in Chapter 4.7, the full calculation in Appendix B. 


7.2.2 Characterization Factor 

The characterization factor can be considered as a way to convert NO x or AIC output into 
equivalent CO 2 mass. The equations are given in Chapter 3.5 and are based around sustained 
global temperature potential and an altitude-dependent forcing factor. The SGTP values can 
be found in Table 3.3 and the forcing factor values in Figure 3.4 or Appendix C. 


7.3 Results Environmental Impact 

7.3.1 Resource Depletion 

Resource depletion is considered equal to fuel consumption as stated before. Fuel consump¬ 
tion results can be found in Chapter 5.3.3. 


7.3.2 Emission Index NO x 

The emission index of NO x is displayed as g of NO x emitted per kg of fuel burned. Calculated 
using BFFM2 (see Chapter 3.5 and Appendix B), the emission index is mostly based on the 
fuel flow of the aircraft. 
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Mach number 


0,4 

0,45 

0,5 

0,55 

0,6 

0,65 

0,7 

0,75 

0,8 


3000 

13,38 

13,63 

14,49 

15,85 

17,72 

20,31 

23,65 

28,19 

35,7 J 


3500 

13,66 

13,75 

14,45 

15,65 

17,35 

19,72 

22,80 

27,00 

33,95 


4000 

14,00 

13,92 

14,46 

15,51 

17,04 

19,21 

22,05 

25,92 

32,35 


4500 

14,39 

14,14 

14,53 

15,42 

16,79 

18,77 

21,38 

24,95 

30,89 


5000 

14,85 

14,42 

14,66 

15,40 

16,60 

18,40 

20,80 

24,09 

29,57 


5500 

15,39 

14,77 

14,84 

15,43 

16,49 

18,11 

20,30 

23,32 

28,39 


6000 

15,99 

15,19 

15,10 

15,54 

16,44 

17,89 

19,88 

22,66 

27,34 


6500 

16,68 

15,68 

15,42 

15,71 

16,46 

17,75 

19,56 

22,10 

26,41 

E 

7000 

17,46 

16,25 

15,82 

15,95 

16,56 

17,69 

19,32 

21,64 

25,60 

CU 

7500 

18,34 

16,90 

16,29 

16,27 

16,73 

17,71 

19,17 

21,28 

24,92 

13 

■M 

8000 

19,33 

17,65 

16,86 

16,68 

16,99 

17,82 

19,11 

21,02 

24,36 

< 

8500 

20,44 

18,51 

17,51 

17,17 

17,34 

18,01 

19,14 

20,87 

23,93 


9000 

21,68 

19,48 

18,27 

17,76 

17,78 

18,30 

19,27 

20,82 

23,61 


9500 

23,06 

20,57 

19,15 

18,46 

18,32 

18,69 

19,51 

20,88 

23,42 


10000 

24,60 

21,80 

20,14 

19,27 

18,97 

19,19 

19,86 

21,06 

23,37 


10500 

26,33 

23,18 

21,27 

20,20 

19,74 

19,81 

20,33 

21,37 

23,44 


11000 

26,26 

23,12 

21,21 

20,14 

19,68 

19,75 

20,27 

21,31 

23,38 

11500 

29,01 

25,39 

23,16 

21,85 

21,20 

21,13 

21,53 

22,46 

24,41 


12000 

32,12 

27,98 

25,37 

23,80 

22,96 

22,73 

23,00 

23,82 

25,64 


12500 

35,65 

30,91 

27,89 

26,02 

24,96 

24,57 

24,70 

25,40 

27,10 


Figure 7.2: Results emission index of NOx in g per kg of fuel 


The general pattern of the EI_NO x results is found in the uncorrected fuel flow results in Fig¬ 
ure 7.3. There is a slight change in some of higher values because of the addition of a humidi¬ 
ty factor for EI_NO x . For a minimal emission of nitrogen oxides, flying lower and slower is 
the best approach. The shape somewhat follows the general shape of the fuel consumption re¬ 
sults but shifted towards lower speed and altitude. 
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0,48 
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0,53 
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3500 
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0,59 

0,65 

0,75 

0,87 

1,04 
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0,59 

0,65 
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0,86 
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1,28 
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0,58 

0,56 

0,57 

0,60 

0,65 

0,73 

0,84 

0,99 

1,24 


5500 

0,61 

0,58 

0,58 

0,61 

0,66 

0,73 

0,83 

0,97 

1,20 


6000 

0,65 

0,61 

0,61 

0,63 

0,67 
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0,83 

0,95 

1,17 


6500 

0,69 
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0,68 

0,74 
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0,94 

1,14 

£ 

7000 

0,74 
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0,75 

0,82 

0,93 

1,12 

CU 

"D 

7500 

0,79 
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0,69 
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0,76 

0,83 

0,93 

1,10 

13 

8000 

0,85 
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0,73 
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0,74 

0,78 
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0,93 

1,09 

’+-> 

< 

8500 

0,92 

0,82 
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0,77 

0,80 

0,85 

0,94 

1,09 


9000 

0,99 

0,88 

0,82 

0,80 
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0,95 

1,09 


9500 

1,08 

0,95 

0,88 

0,85 

0,84 
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0,90 

0,97 

1,10 


10000 

1,18 

1,03 

0,95 

0,90 

0,89 

0,90 

0,93 

1,00 

1,11 


10500 

1,29 

1,12 

1,02 

0,97 

0,94 

0,95 

0,97 

1,03 

1,14 


11000 

1,30 

1,14 

1,04 

0,98 

0,96 

0,96 

0,99 

1,04 

1,15 


11500 

1,51 

1,31 

1,19 

1,12 

1,08 

1,08 

1,10 

1,15 

1,26 


12000 

1,76 

1,52 

1,37 

1,28 

1,23 

1,22 

1,23 

1,28 

1,39 

12500 

2,05 

1,76 

1,58 

1,47 

1,41 

1,38 

1,39 

1,43 

1,53 


Figure 7.3: Results corrected fuel flow in kg per second 
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7.3.3 Equivalent CO 2 Mass 


The equivalent CO 2 mass is represented as kg of CO 2 emitted per NM flown, per seat. By di¬ 
viding by the number of seats, the metric becomes more generic. 
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0,109 
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0,146 
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a; 

-Q 
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0,184 

0,200 

13 

■M 

8000 
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0,236 

0,252 

-t— > 

< 

8500 

0,458 

0,391 
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0,310 

0,315 

0,332 


9000 

0,587 

0,502 

0,453 

0,424 

0,407 
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0,397 
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0,419 


9500 

0,710 

0,607 

0,547 

0,512 

0,490 

0,479 

0,475 

0,478 

0,494 


10000 

0,852 

0,723 

0,649 

0,605 

0,578 

0,563 

0,556 

0,557 

0,572 


10500 

0,975 

0,806 

0,710 

0,653 

0,618 

0,597 

0,587 

0,586 

0,601 


11000 

1,003 

0,810 

0,702 

0,638 

0,599 

0,576 

0,565 

0,565 

0,582 


11500 

1,123 

0,873 

0,733 

0,650 

0,599 

0,569 

0,553 

0,550 

0,566 


12000 

1,270 

0,939 

0,755 

0,646 

0,579 

0,539 

0,516 

0,509 

0,524 


12500 

1,504 

1,079 

0,843 

0,704 

0,618 

0,566 

0,536 

0,523 

0,535 


Figure 7.4: Results equivalent C02 mass in kg per NM flown, per seat 


There is a relatively clear color transition around 7 000 m (about 23 000 ft) in these results. 
This is the point at which the influence of aircraft-induced clouds starts to rise rapidly and 
where the full warming effect of NO x becomes available. The influence of AIC appears in the 
form of the forcing factor s and can be seen on Figure 3.4. On Figure 7.4, a vaguely visible 
orange colored band can be seen around 10 000 to 11 500 m (about 33 000 to 38 000 ft). This 
is the altitude associated with a maximum value of the forcing factor for AIC, and thus the 
biggest contribution. It also happens to be the most common cruise altitude zone for passen¬ 
ger jet aircraft nowadays. 


Furthermore, in the green zone at lower altitude and moderate speed a shape resembling the 
fuel consumption pattern can be seen. This is because fuel consumption directly influences 
the amount of CO 2 and NO x emitted (see Eq. (3.17)). 


Equivalent CO 2 mass (Eq. (3.17)) is a rather complex metric that depends on numerous de¬ 
pendent and independent parameters. To get a better understanding of the influence of differ¬ 
ent emission contributors on the equivalent CO 2 mass, a comparison of different flight param¬ 
eter scenarios is added below. 
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Figure 7.5: Importance of individual parameters in the equivalent C02 mass equation 


The first observation is that aircraft-induced clouds have zero or close to zero percent influ¬ 
ence in the lower altitude zone. The reason being the low forcing factor values for AIC (see 
Figure 3.4). At higher altitudes the influence of AIC contributes for the majority of the envi¬ 
ronmental impact. 


Another important note is that NOx has a negative contribution to the equivalent CO 2 mass at 
low altitudes. Nitrogen oxides can cause both a warming and a cooling effect. The warming 
effect occurs because of short ozone (O3) and the cooling effect because of long ozone and 
methane (O3 and CH4 respectively). This phenomenon is visible in a cause effect fonn in Ta¬ 
ble 2.1 and in a quantitative form as SGTP values in Table 3.3. This cooling effect (or nega¬ 
tive contribution) stops at around 6 000 m (about 19 500 ft). 


7.3.4 Environmental Impact 

The environmental impact is a combination of resource depletion and engine emissions. The 
impact of engine emissions is made quantitative by use of the equivalent CO 2 mass metric in 
the previous subchapter. Both resource depletion and engine emissions are normalized, multi- 
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plied by their respective weighting factor and added up to form the environmental impact. 
Weighting factors can be altered to change the importance of each contributor. Altering the 
weighting factors will result in a different outcome. The value of these weighting factors 
could be a topic of discussion, but for this result a neutral approach with 50 - 50 weighting 
factors (and thus an equal importance of resource depletion and engine emissions) is used. 
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Figure 7.6: Results environmental impact (normalized between 0 and 1) 



□ 0 , 000 - 0,100 □ 0 , 100 - 0,200 □ 0 , 200 - 0,300 □ 0 , 300 - 0,400 □ 0 , 400 - 0,500 

□ 0 , 500 - 0,600 □ 0 , 600 - 0,700 □ 0 , 700 - 0,800 □ 0 , 800 - 0,900 10 , 900 - 1,000 


Figure 7.7: 


Contour plot of the environmental impact results 
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Due to these results being a combination of fuel consumption and equivalent CO 2 mass, a lot 
of similarities are found with both models. First, a clear altitude boundary can be seen where 
the green zone turns into yellow. The altitude has shifted from around 7 000 meter in the orig¬ 
inal equivalent CO 2 mass model to 8 000 meter here due to the shape of the fuel consumption 
results. Furthermore, the diagonal shape of the fuel consumption pattern from Chapter 5.3.3 is 
visible. This implies a reduced environmental impact at higher altitudes when flying faster, 
although this effect is rather limited. In conclusion, the results come down to flying signifi¬ 
cantly lower and slightly slower for a minimal environmental impact. 


7.4 Conclusion Environmental Impact 

In this chapter the environmental impact of an aircraft is represented in a quantitative way by 
combining resource depletion and engine emissions. Resource depletion is calculated as fuel 
consumption and engine emissions by way of equivalent CO 2 mass. The equivalent CO 2 mass 
is dependent on varying flight parameters by means of emission indices and altitude- 
dependent characterization factors. It combines the effects from CO 2 and non-CCfi emissions 
into a single metric. 

There exists a certain level of uncertainty for several contributors to the environmental impact 
of aircraft. Especially the influence of aircraft-induced cloudiness (AIC) in the form of con¬ 
trails and cirrus clouds has a poor level of scientific understanding (see Figure 2.5). This and 
the fact that seasonal and location-based influences on contrail fonnation have been neglect¬ 
ed, results in some inaccuracy in the final environmental impact results. This does not pro¬ 
duce any issues when solely comparing results with each other to achieve a trend in environ¬ 
mental impact variation. 

The main conclusion of the results in Figure 7.6 is that a reduction in cruise altitude would 
greatly reduce the environmental impact of aircraft. Looking at the equivalent CO 2 mass re¬ 
sults (Figure 7.4), the common cruise altitude zone for jet aircraft nowadays (10 500 m to 
11 500 m) results in the biggest equivalent CO 2 mass emission. Flying lower would greatly 
reduce the fonnation of contrails. In Figure 3.4 can be seen that the forcing factor of AIC 
shifts towards zero below 6 500 m, meaning no more fonnation of contrails and aircraft- 
induced cirrus clouds. 

The problem with a lower cruise altitude is the decrease in overall engine efficiency and the 
corresponding increase in TSFC. Furthermore, the aircraft’s aerodynamic efficiency is not op¬ 
timized for this lower altitude. Meaning that it would be far out of the optimal dark green 
zone of the aerodynamic efficiency results in Figure 5.5, which in turn increases fuel con¬ 
sumption. A smaller wing surface resulting from an increase in wing loading will optimize the 
aircraft’s aerodynamic efficiency for lower cruise altitudes. This scenario is analyzed in Chap¬ 
ter 8.3. 
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From an operational standpoint flying at lower altitudes also proves some challenges. A first 
issue relates to the historical reason of flying at higher altitudes: when flying at or around 
tropopause level an aircraft can fly over most weather phenomena, reducing turbulence and 
flying through thunderstonns. This is both a safety and passenger comfort measure. Addition¬ 
ally, bird strikes are extremely rare at tropopause level (90 % occurs under 3 000 ft according 
to Cleary 2005), further increasing safety at higher altitudes. Another reason is the added 
safety margin for an engine failure or other aircraft defect at high altitude. The flight crew has 
more time to follow necessary procedures, communicate, or come up with other measures if 
needed. 

Changing the regular cruise altitude of an Airbus A320-200 of about 11 500 m to an altitude 
of 6 500 m at a constant Mach 0.78 would result in a decrease of equivalent CO 2 mass of 
77.66 % and an increase of fuel consumption of 5.64 %. The increase of fuel consumption is 
mostly influenced by an increase of TSFC of 6.04 % and a decrease of the aerodynamic effi¬ 
ciency of 5.44 %. Combining equivalent CO 2 mass and resource depletion (fuel consumption) 
into the environmental impact would result in a decrease of 69.59 % in environmental impact. 

A secondary conclusion is that the cruise Mach number must be reduced in order to minimize 
the environmental impact. Flying slower does not have a substantial impact on the equivalent 
CO 2 mass alone, but when combining the need of flying lower with resource depletion (fuel 
consumption results), it is clear that a reduction in cruise altitude requires a reduction in cruise 
speed in order to maintain minimum fuel consumption. 

At an altitude of 6 500 m a cruise speed reduction to Mach 0.65 would result in a minimal en¬ 
vironmental impact. Resulting in a further reduction of 82.61 % in environmental impact 
compared to flying at Mach 0.78 at the same altitude. This can be explained by the fuel con¬ 
sumption pattern in Figure 5.10: Mach 0.65 at 6 500 m is right in the optimal dark green zone 
while flying Mach 0.78 at this altitude results in an increase in fuel consumption of 11.86 %. 
The emitted equivalent CO 2 mass is reduced by 13.36 % with this cruise speed reduction. 

The altitude of 6 500 m is used as an example in the previous paragraphs because of the lack 
of contrail fonnation and the fact that this is the optimum altitude for the combined model in 
Chapter 8. The actual flight parameters for a minimal environmental impact are Mach 0.6 at 
an altitude of 5000 m. Comparing this to the regular values of Mach 0.78 at 11 500 m, a re¬ 
duction in environmental impact of 98.81 % would be possible. This is however a highly un¬ 
realistic flight scenario. The lower altitude poses problems as stated before and an even lower 
cruise Mach number results in a greater economic loss. 

There is a range of other options to decrease the environmental impact of operating an aircraft 
apart from changing its flight parameters. New engine design to lower the thrust specific fuel 
consumption and thus the emission of CO 2 and NO x is an option. Small improvements like 
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adding an electric motor to the nose wheel for ground handling at airports instead of using the 
aircraft engines for this (DLR 2011). Using biofuel or a mixture of biofuel and regular jet fuel 
is another option which is already being tested by multiple manufacturers, among them Boe¬ 
ing with their ecoDemonstrator program (Jensen 2018). Even fully electric or hybrid electric 
aircraft designs have been considered but will not be viable in the near future (Scholz 2018). 
The drawback of all these options is that a modification to the existing aircraft or a new air¬ 
craft design is necessary, or that they are simply not viable yet (e.g.: biofuel or electric air¬ 
craft). Having to invest in a new aircraft or extensive modifications makes all these options 
less attractive to airlines due to the long lifespan of their current aircraft. In contrast, small 
changes to flight parameters can be implemented right now without major challenges. 
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8 Flight Parameter Optimization: Combined Model 

8.1 General 

The combined model is a combination of seat-mile cost to represent the economic conse¬ 
quences and environmental impact for ecologic consequences. Fuel consumption will not be 
used as a direct input in the model because it is already heavily implemented in both DOC and 
environmental impact models. 

Fuel consumption is present in DOC in the fonn of the Breguet factor, which determines the 
necessary fuel mass needed to complete the flight. It influences the environmental impact 
model in three ways. First, directly in the fonn of resource depletion, responsible for half of 
the environmental impact results (when using the default 50 - 50 weighting factors). Second, 
it is also used to calculate the aircraft’s fuel flow, necessary for calculating the emission index 
of NOx. Finally, fuel consumption is used to calculate the equivalent CO 2 mass with Eq. 
(3.17). 



Figure 8.1: Flowchart combined model calculation 


8.2 Results Combined Model 

The results for the combined model can vary widely depending on the distribution of 
weighting factors. Two scenarios are presented here: a perhaps idealistic approach with a 50 - 
50 weight distribution between economic and ecologic importance and a perhaps more realis- 
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tic approach of 80 % economic importance and 20 % ecologic importance. These are only ar¬ 
bitrary values and can be changed in the Excel tool to generate different results. 


8.2.1 Equal Importance DOC and Environmental Impact Case 


In this case the combined model is a result with equal weighting factors of 50 % for DOC and 
50 % for environmental impact. 



Mach number 


0,4 

0,45 

0,5 

0,55 

0,6 

0,65 

0,7 

0,75 

0,8 

3000 

0,381 

0,271 

0,196 

0,141 

0,101 

0,075 

0,060 

0,056 

0,078 

3500 

0,401 

0,287 

0,208 

0,151 

0,109 

0,080 

0,062 

0,056 

0,073 

4000 

0,413 

0,295 

0,213 

0,154 

0,110 

0,080 

0,060 

0,051 

0,065 

4500 

0,426 

0,303 

0,219 

0,157 

0,111 

0,079 

0,057 

0,047 

0,057 

5000 

0,441 

0,314 

0,225 

0,161 

0,114 

0,079 

0,056 

0,043 

0,050 

5500 

0,459 

0,326 

0,234 

0,168 

0,118 

0,081 

0,056 

0,041 

0,045 

6000 

0,477 

0,339 

0,243 

0,173 

0,121 

0,083 

0,056 

0,039 

0,040 

6500 

0,498 

0,354 

0,253 

0,181 

0,126 

0,086 

0,056 

0,038 

0,036 

E " 7000 

0,527 

0,376 

0,270 

0,194 

0,137 

0,094 

0,063 

0,042 

0,038 

7 i 7500 

“0 

0,558 

0,400 

0,288 

0,209 

0,149 

0,104 

0,071 

0,048 

0,041 

2 8000 

0,595 

0,429 

0,312 

0,229 

0,166 

0,118 

0,082 

0,057 

0,048 

3 = 8500 

0,635 

0,466 

0,342 

0,254 

0,187 

0,137 

0,098 

0,071 

0,060 

9000 

0,680 

0,507 

0,376 

0,282 

0,212 

0,158 

0,117 

0,088 

0,074 

9500 

0,726 

0,548 

0,410 

0,315 

0,243 

0,189 

0,148 

0,117 

0,103 

10000 

0,777 

0,592 

0,447 

0,346 

0,270 

0,213 

0,169 

0,136 

0,119 

10500 

0,828 

0,629 

0,480 

0,372 

0,291 

0,230 

0,183 

0,148 

0,128 

11000 

0,839 

0,635 

0,483 

0,373 

0,291 

0,229 

0,182 

0,147 

0,128 

11500 

0,882 

0,662 

0,505 

0,388 

0,301 

0,236 

0,186 

0,149 

0,127 

12000 

0,932 

0,691 

0,526 

0,401 

0,309 

0,239 

0,186 

0,146 

0,122 

12500 

1,000 

0,734 

0,556 

0,428 

0,328 

0,254 

0,197 

0,154 

0,128 

Figure 8.2: 

Results combined 

model with equal importance of DOC and environmental 

impact 


(normalized between 0 and 1) 
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□ 0 , 000 - 0,100 □ 0 , 100 - 0,200 □ 0 , 200 - 0,300 □ 0 , 300 - 0,400 □ 0 , 400 - 0,500 

□ 0 , 500 - 0,600 □ 0 , 600 - 0,700 □ 0 , 700 - 0,800 □ 0 , 800 - 0,900 □ 0 , 900 - 1,000 


Figure 8.3: Contour plot of the combined model results (DOC and environmental impact) 

The influence of the DOC is clear in the way that it shifts the optimum for minimum envi¬ 
ronmental impact towards the right, to a higher cruise Mach number. The shape from Figure 
7.6 of the environmental impact results is still clearly visible with its tendency towards lower 
altitude. These two combine to an optimum result at a lower altitude and an unchanged cruise 
Mach number compared to regular flight parameters. 
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8.2.2 Emphasis on DOC Case 


In this case the combined model is a result with a weighting factor of 80 % for DOC and 20 % 
for environmental impact. 


Mach number 



0,4 

0,45 

0,5 

0,55 

0,6 

0,65 

0,7 

0,75 

0,8 


3000 

0,577 

0,420 

0,306 

0,219 

0,151 

0,099 

0,061 

0,035 

0,031 


3500 

0,604 

0,442 

0,326 

0,236 

0,166 

0,112 

0,071 

0,043 

0,036 


4000 

0,618 

0,452 

0,333 

0,242 

0,171 

0,116 

0,073 

0,043 

0,034 


4500 

0,632 

0,463 

0,341 

0,248 

0,175 

0,119 

0,075 

0,043 

0,031 


5000 

0,648 

0,474 

0,350 

0,255 

0,181 

0,122 

0,077 

0,044 

0,029 


5500 

0,666 

0,488 

0,359 

0,263 

0,186 

0,127 

0,080 

0,045 

0,027 


6000 

0,684 

0,502 

0,369 

0,270 

0,192 

0,130 

0,082 

0,045 

0,025 


6500 

0,704 

0,517 

0,380 

0,279 

0,199 

0,135 

0,085 

0,046 

0,024 

E 

7000 

0,728 

0,535 

0,394 

0,290 

0,207 

0,142 

0,089 

0,049 

0,025 

(U 

7500 

0,754 

0,555 

0,409 

0,302 

0,217 

0,149 

0,095 

0,053 

0,026 

=5 

+-> 

8000 

0,783 

0,579 

0,428 

0,316 

0,229 

0,159 

0,102 

0,058 

0,029 

< 

8500 

0,807 

0,605 

0,449 

0,333 

0,242 

0,169 

0,111 

0,065 

0,034 


9000 

0,833 

0,635 

0,472 

0,353 

0,258 

0,182 

0,122 

0,073 

0,040 


9500 

0,860 

0,665 

0,497 

0,378 

0,285 

0,211 

0,151 

0,103 

0,070 


10000 

0,890 

0,694 

0,523 

0,401 

0,304 

0,227 

0,165 

0,115 

0,080 

10500 

0,919 

0,717 

0,549 

0,421 

0,321 

0,241 

0,177 

0,125 

0,088 


11000 

0,931 

0,726 

0,557 

0,426 

0,325 

0,245 

0,180 

0,128 

0,091 


11500 

0,950 

0,738 

0,572 

0,438 

0,334 

0,251 

0,184 

0,131 

0,092 


12000 

0,971 

0,750 

0,586 

0,449 

0,341 

0,255 

0,187 

0,131 

0,090 

12500 

1,000 

0,769 

0,599 

0,466 

0,354 

0,266 

0,195 

0,137 

0,095 


Figure 8.4: Results combined model with an emphasis on DOC (normalized between 0 and 1) 



□ 0 , 000 - 0,100 □ 0 , 100 - 0,200 □ 0 , 200 - 0,300 □ 0 , 300 - 0,400 □ 0 , 400 - 0,500 

□ 0 , 500 - 0,600 □ 0 , 600 - 0,700 □ 0 , 7000,800 □ 0 , 800 - 0,900 □ 0 , 900 - 1,000 


Figure 8.5: Contour plot of the combined model results (DOC) 
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In these results the shape of the environmental impact results is no longer clearly visible. All 
focus is aimed towards an economic optimum, which is achieved at the highest Mach number. 
There is still a slight shift towards the optimal altitude for environmental impact, but less visi¬ 
ble than in the previous set of results in Chapter 8.2.1. 

Another observation is that the normalized values around the optimum zone are lower than in 
the previous case. This is because the influence of environmental impact shifts the optimum 
zone towards a smaller Mach number. And with this influence less present than before, the 
optimum zone is mostly dictated by the DOC. Optimum DOC is at a maximum Mach num¬ 
ber, therefore the values come closer to an optimum normalized value of zero. 


8.3 Influence of the Wing Loading 

In this chapter the combined model from Chapter 8.2 will be altered with a higher wing load¬ 
ing. Wing loading is defined as the ratio of an aircraft’s mass (MTOW) over its wing surface. 
Wing loading and thrust-to-weight ratio are the two most important parameters to be opti¬ 
mized during preliminary aircraft design. Both are represented graphically in a matching chart 
with the aim to achieve the highest possible wing loading and lowest possible thrust-to-weight 
ratio. As can be seen in Figure 8.6, maximum wing loading is limited by the available landing 
field length and requires an increasing thrust-to-weight ratio for take-off (Scholz 2015). 



Figure 8.6: Illustrative example of a matching chart (Scholz 2015) 
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Wing loading influences the lift coefficient, as given in the following equation where the lift 
coefficient is represented in function of wing loading: 


C L (m/5) = 


2 • g • (m/S) 
p ■ M 2 ■ a 2 


( 8 . 1 ) 


This equation reveals that for a fixed altitude (represented by p and a) a higher wind loading 
will produce a higher lift coefficient. This optimizes the aerodynamic efficiency for lower al¬ 
titudes, therefore shifting the zone of minimal fuel consumption to a lower altitude. The main 
purpose of choosing a higher wing loading for an aircraft is to lower its wing surface for the 
same MTOW. A smaller wing surface is optimal for flying at lower altitudes. 

An increase in wing loading also means a required increase in landing field length. Landing 
field length is defined as the landing distance multiplied by a safety factor (1.667 for jet air¬ 
craft), as can be seen in Figure 8.7. 



< - Landing Distance -H 

4 - Landing Field Length -► 

Figure 8.7: Landing field length (Scholz 2015) 

The landing field length s LFL relates to the ratio of maximum landing mass over wing surface 
in the following way: 


m ML 

~S~ 


— k, - a ■ C, 


L , max,L " ^LFL 


( 8 . 2 ) 


From this ratio the wind loading can be calculated by means of the maximum landing mass to 
maximum take-off mass ratio: 


m MTOW _ m ML /$ ^ 2 ) 

5 m ML / m MTO 

When assuming the maximum lift coefficient during landing as constant, the landing field 
length is directly proportional with wing loading. 
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As an experiment, the original wing loading of the Airbus A320-200 of 600.49 kg/m 2 is in¬ 
creased with about 50 % to 900.74 kg/m 2 . As a result of Eq. (8.2) and (8.3), the landing field 
length increases with 50 % as well. This makes it impossible for this aircraft to land at many 
airports where it normally would have been able to. 


This 50 % increase in wing loading yields the following combined model results at 50 - 50 
weighting factors for DOC and environmental impact: 


Mach number 
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0,45 

0,5 
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3500 
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0,308 

0,225 
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0,114 

0,077 

0,049 

0,030 

0,024 


4000 

0,441 

0,319 

0,233 

0,169 

0,119 

0,081 

0,052 

0,032 

0,024 


4500 

0,457 

0,331 

0,241 

0,175 

0,124 

0,085 

0,055 

0,033 

0,024 


5000 

0,474 

0,343 

0,251 

0,183 

0,131 

0,090 

0,059 

0,036 

0,024 


5500 

0,492 

0,358 

0,262 

0,192 

0,138 

0,096 

0,063 

0,039 

0,026 


6000 

0,507 

0,373 

0,274 

0,201 

0,145 

0,101 

0,067 

0,041 

0,027 


6500 

0,524 

0,391 

0,287 

0,211 

0,153 

0,107 

0,072 

0,045 

0,029 

E 
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0,546 

0,412 

0,303 

0,225 
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0,117 

0,079 

0,051 

0,033 

(U 
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0,569 

0,433 

0,322 
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0,177 
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0,088 

0,058 

0,039 

=3 

■M 
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0,596 

0,455 

0,344 

0,258 

0,192 

0,140 
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0,067 

0,046 

< 

8500 

0,630 

0,481 

0,370 

0,280 

0,210 

0,156 

0,113 

0,079 

0,056 


9000 

0,668 

0,511 

0,398 

0,304 

0,231 

0,174 

0,129 

0,093 

0,068 


9500 

0,707 

0,541 

0,423 

0,333 

0,259 

0,201 

0,156 

0,120 

0,094 


10000 

0,752 

0,576 

0,450 

0,362 

0,283 

0,223 

0,174 

0,137 

0,109 


10500 

0,802 

0,611 

0,477 

0,384 

0,305 

0,241 

0,191 

0,151 

0,121 


11000 

0,817 

0,621 

0,484 

0,388 

0,308 

0,243 

0,192 

0,151 

0,122 


11500 

0,866 

0,653 

0,506 

0,404 

0,325 

0,256 

0,202 

0,160 

0,129 


12000 

0,926 

0,691 

0,531 

0,421 

0,340 

0,270 

0,212 

0,167 

0,134 


12500 

1,000 

0,739 

0,565 

0,446 

0,359 

0,290 

0,230 

0,182 

0,147 


Figure 8.8: Results combined model for a 50 % higher wing loading 
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Figure 8.9: Contour plot of the combined model results (50 % higher wing loading) 
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In the following figures, the influence of a higher wing loading on fuel consumption, DOC, 
and environmental impact will be presented. 



Mach number 


0,4 

0,45 

0,5 

0,55 

0,6 

0,65 

0,7 

0,75 

0,8 


3000 

9,76E-05 

8,14E-05 

7,24E-05 

6,76E-05 

6,57E-05 

6,64E-05 

6,89E-05 

7,39E-05 

8,45E-05 


3500 

l,02E-04 

8,44E-05 

7,42E-05 

6,86E-05 

6,59E-05 

6,59E-05 

6,78E-05 

7,20E-05 

8,15E-05 


4000 

l,08E-04 

8,79E-05 

7,64E-05 

6,98E-05 

6,64E-05 

6,57E-05 

6,70E-05 

7,05E-05 

7,90E-05 


4500 

l,13E-04 

9,18E-05 

7,90E-05 

7,14E-05 

6,73E-05 

6,59E-05 

6,64E-05 

6,92E-05 

7,68E-05 


5000 

l,20E-04 

9,62E-05 

8,19E-05 

7,34E-05 

6,84E-05 

6,63E-05 

6,62E-05 

6,84E-05 

7,49E-05 


5500 

l,27E-04 

l,01E-04 

8,54E-05 

7,57E-05 

6,99E-05 

6,71E-05 

6,63E-05 

6,78E-05 

7,35E-05 


6000 

l,35E-04 

l,07E-04 

8,93E-05 

7,84E-05 

7,17E-05 

6,82E-05 

6,68E-05 

6,76E-05 

7,24E-05 


6500 

l,44E-04 

l,13E-04 

9,37E-05 

8,16E-05 

7,40E-05 

6,96E-05 

6,75E-05 

6,77E-05 

7,17E-05 


7000 

l,54E-04 

l,20E-04 

9,87E-05 

8,52E-05 

7,66E-05 

7,15E-05 

6,87E-05 

6,81E-05 

7,13E-05 

cu 

-0 

7500 

l,65E-04 

l,27E-04 

l,04E-04 

8,94E-05 

7,97E-05 

7,37E-05 

7,02E-05 

6,90E-05 

7,14E-05 

13 

■M 

8000 

l,77E-04 

l,36E-04 

l,llE-04 

9,42E-05 

8,33E-05 

7,64E-05 

7,21E-05 

7,02E-05 

7,18E-05 

4— 1 

< 

8500 

l,90E-04 

l,45E-04 

l,18E-04 

9,95E-05 

8,74E-05 

7,95E-05 

7,45E-05 

7,18E-05 

7,26E-05 


9000 

2,04E-04 

l,56E-04 

l,25E-04 

l,06E-04 

9,21E-05 

8,32E-05 

7,73E-05 

7,39E-05 

7,39E-05 


9500 

2,21E-04 

l,68E-04 

l,34E-04 

l,12E-04 

9,75E-05 

8,74E-05 

8,06E-05 

7,65E-05 

7,56E-05 


10000 

2,39E-04 

l,81E-04 

l,44E-04 

l,20E-04 

l,03E-04 

9,22E-05 

8,45E-05 

7,95E-05 

7,79E-05 


10500 

2,59E-04 

l,95E-04 

l,55E-04 

l,28E-04 

l,10E-04 

9,77E-05 

8,90E-05 

8,31E-05 

8,06E-05 


11000 

2,58E-04 
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Figure 8.10: Results fuel consumption in kg fuel per NM flown, per kg of aircraft MTOW (50% 
higher wind loading) 


The smaller wing surface that results from a higher wing loading will shift this aircraft’s aero¬ 
dynamic efficiency optimum to a lower altitude. Which consequently shifts the minimal fuel 
consumption optimum to a lower altitude. There will be a small increase in TSFC at lower al¬ 
titude, but the improved aerodynamic efficiency outweighs this for the fuel consumption re¬ 
sult. 
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Figure 8.11: Results seat-mile cost in USD (50 % higher wing loading) 


The seat-mile cost experiences a slight change by increasing wing loading. Seat-mile cost at 
lower altitudes decreases while it increases at higher altitudes compared to the regular wing 
loading case. Seat-mile cost depends on the annual DOC of an aircraft. The annual DOC in 
turn depends on fuel consumption results and the optimal zone here is shifted towards lower 
altitudes. 
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Figure 8.12: Results environmental impact per NM flown per seat (50 % higher wing loading) 


Comparing to the environmental impact with a regular wing loading, a decrease of the opti¬ 
mum value is visible. The optimal value with regular wing loading is 0.0022 normalized be¬ 
tween 0 and 1. The optimal value for a higher wing loading is 3.74 ■ 10~ 6 , also normalized 
between 0 and 1. This value is significantly smaller because the zones of minimal fuel con¬ 
sumption and minimal environmental impact start to merge with a higher wing loading. Es¬ 
sentially, the higher wing loading brings fuel consumption and equivalent CO 2 mass optima 
closer to each other. 


8.4 Conclusion Combined Model 

In this chapter all three separate models have been combined into one model. The goal is to 
find optimum flight parameters for both economic and ecologic interests. The combined mod¬ 
el is calculated for two cases: an equal importance of economic and ecologic benefits and a 
case with emphasis on the economic benefits. 

The accuracy of this model is only as strong as the combined accuracies of all separate mod¬ 
els. As mentioned before, all individual results are calculated with the same methods. There¬ 
fore, a comparison of results is still accurate enough to visualize a trend with changing flight 
parameters. 

The general trend in both cases is a reduction in cruise altitude. This has a positive effect on 
seat-mile cost because of the increased TAS and on the environmental impact because of the 
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reduced equivalent CO 2 mass. A decrease in cruise Mach number on the other hand would 
benefit the environmental impact but would impose a penalty on the seat-mile cost. An im¬ 
portant point to consider is that flying fast and low could be outside of the aircraft’s permissi¬ 
ble flight envelope (see Chapter 9.2). 

In the 50 - 50 DOC and environmental model, the optimum combined value can be found at 
Mach 0.8 at 6 500 m. This results in a combined optimization of 71.30 % compared to the 
regular flight parameters (Mach 0.78 at 11 500 m). Mach 0.8 is close to the maximum operat¬ 
ing Mach number of the aircraft but results in the optimum seat-mile cost. Corresponding to 
this speed, an altitude of 6 500 m results in the lowest environmental impact (see Figure 7.6). 
In the case of an emphasis on DOC, the optimum value is still at the same flight parameters 
but has an ever lower normalized value. This is because the optimal DOC zone is already at 
this point. 

A higher wing loading implies a smaller wing surface when keeping the MTOW constant. 
This results in a shift of the minimum fuel consumption zone towards lower altitude, moving 
it closer to the zone of minimal equivalent CO 2 mass. An increase in wing loading will essen¬ 
tially merge the zones of minimum fuel consumption and minimum equivalent CO 2 mass. As 
a result, the normalized environmental impact values are significantly lower than in the nor¬ 
mal wing loading case. 

The optimal combined value with a 50 % higher wing loading can be found at Mach 0.8 at 
3 000 m. This result in a combined optimization of 84.32 %, an increase of 13.02 % compared 
to the normal wing loading case. This is however at a greatly reduced altitude and will proba¬ 
bly be outside of the aircraft’s pennissible flight envelope as mentioned before. Besides the 
flight envelope problem, flying at low altitudes generates safety and passenger comfort issues 
(see Chapter 7.4). 

An important recent development that could alter the relationship between DOC and envi¬ 
ronmental impact is a CO 2 emissions standard for aircraft manufacturers. This standard is 
drafted by the ICAO and will be implemented starting in 2020 (scheduled date at the moment 
of publication). In short, the standard applies to all aircraft with an MTOW over 5 700 kilo¬ 
gram, which implies most aircraft types in the current global passenger jet fleet. It is based on 
fuel consumption during the cruise flight stage only. The standard covers three categories: 
newly designed aircraft which require a type certificate will have to comply by January 2020, 
all aircraft in production which require modifications will have to comply by the end of 2023, 
and finally all remaining in-production aircraft will have to comply by 2028 (ATAG 2016). 
The implication of this standard could be a shift of importance from economic to ecologic 


reasons. 
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9 Conclusions and Recommendations 

9.1 Conclusions 

The main goal of this work is to visualize and analyze the influence of cruise Mach number, 
cruise altitude, and aircraft wing loading on fuel consumption, DOC, and environmental im¬ 
pact. And subsequently combine these three objectives into one model. All calculations are 
done on the Airbus A320-200 as a reference aircraft, but the models are intended for any pas¬ 
senger jet aircraft. 

Fuel consumption - Minimum fuel consumption is achieved by following a rising Mach 
number and rising altitude pattern. Each Mach number has its corresponding optimal altitude. 
A low Mach number must be flown at lower altitudes and an increasing Mach number re¬ 
quires an increase in altitude. Fuel consumption is predominantly influenced by the aerody¬ 
namic efficiency of the aircraft and to a lesser extent by the TSFC, Mach number and local 
speed of sound. A higher wing loading shifts the zone for minimal fuel consumption to lower 
altitudes due to a shift of the optimal lift coefficient (and aerodynamic efficiency) in the same 
direction. 

Direct operating cost - Two metrics are available for calculating an aircraft’s DOC: annual 
DOC and seat-mile cost. The annual DOC is minimal at minimal TAS (small Mach number, 
large altitude). However, this result does not represent the economic impact: it simply mini¬ 
mizes the amount of flights per year which reduces all time-dependent costs (fuel, mainte¬ 
nance, staff, and fees and charges), in this way also minimizing possible revenue generation 
by the aircraft. The more correct approach is to compare the seat-mile cost. 

The seat-mile cost must be as low as possible to increase possible revenue gain. It is predomi¬ 
nantly influenced by the aircraft’s flight time. A shorter flight time implies a lower seat-mile 
cost, situating the optimum at the highest Mach number combined with a medium-low alti¬ 
tude. The higher wing loading case shifts optimal seat-mile cost from medium-low altitudes to 
the lowest altitude because of the change in annual DOC through fuel consumption. 

Environmental impact - The environmental impact is expressed as a combination of re¬ 
source depletion (equal to fuel consumption) and engine emissions (in the form of equivalent 
CO 2 mass). The emission of equivalent CO 2 mass peaks at the cruise altitude zone used by 
passenger jet aircraft nowadays. A reduction of cruise altitude from 11 500 m to 6 500 m 
would reduce the fonnation of AIC to almost none. For the Airbus A320-200, this and other 
beneficial effects of flying at a lower altitude (e.g. lower emission of NO x ) would reduce 
equivalent CO 2 mass by 77.66 % with an increase in fuel consumption of 5.64 %. If a varia¬ 
tion in Mach number would also be applied from Mach 0.78 to Mach 0.65 at 6 500 m, a fur- 
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ther reduction of equivalent CO 2 mass of 13.36 % and a reduction in fuel consumption of 
11.86 % are achieved. 

A higher wing loading optimizes the normalized environmental impact results even more. 
This is because the zones of minimum fuel consumption and minimum equivalent CO 2 mass 
merge with increasing wing loading. 

Combined model - In the combined model results are generated for two cases: one where 
economic and ecologic benefits have an equal importance and one with an emphasis on the 
economic benefits. In both cases the general trend is a reduction in cruise altitude and an unal¬ 
tered cruise Mach number compared to present-day conditions. An emphasis on economic in¬ 
terest yields lower (and thus better) results in the optimal zone than the case with equal im¬ 
portance. This is because the optimum is already mostly influence by DOC. Therefore, de¬ 
creasing the importance of environmental impact will further optimize the results. A higher 
wing loading shifts the optimal zone towards lower altitudes. 

Operational limitations - Flying at a lower cruise altitudes comes with certain complica¬ 
tions. As mentioned in Chapter 7.4 most of these are related to aircraft safety or passenger 
comfort: no possibility to fly over weather phenomena, a higher chance of bird strikes, and 
less reaction time in case of emergency. When coupled with low Mach number, there is the 
added possibility that the aircraft is flying outside of its pennitted flight envelope. Load fac¬ 
tors here can be higher than the structural limit of the aircraft or it could be flying below stall 
speed. 


9.2 Recommendations 

When considering varying cruise Mach number and altitude for an aircraft, the flight envelope 
should be kept in mind. In this thesis optimal theoretical values are proposed for several 
goals, but there is no guarantee that the aircraft could actually operate at these flight condi¬ 
tions. As an example: flying at low Mach number and high altitude could be well below the 
aircraft’s stall speed and flying at high Mach number and low altitude could result in ambient 
conditions above the aircraft’s structural limit (the cabin pressure differential would be too 
high). Therefore, a protection could be built into all models that ensures the results to be with¬ 
in the considered aircraft’s flight envelope. 

As mentioned in Chapter 6.2.2, flight time has a considerable influence on the seat-mile cost 
and various DOC categories. This means that an improvement on the approximating flight 
time model introduced in Chapter 6.2.2 would yield an improved accuracy in both seat-mile 
cost and annual DOC results. The model could be enhanced by including BADA data for 
more and more modern aircraft (more up-to-date BADA data is not publicly available at the 
time of publication) and by modelling the data with more fitting and complex equations. 
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Finally, the influence of the ICAO’s CO 2 standard (as described in Chapter 8.4) or possible 
future environmental taxes could be included in a quantitative way in the weighting factors of 
DOC and environmental impact. This standard could shift importance from economic to eco- 
logic consequences or change ecologic into economic consequences. 
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Appendix A - TSFC Model by Herrmann 

All equations are established by Herrmann 2010. 
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T(h) is the temperature at altitude, T 0 = 288 K, T m is the take-off thrust of one engine and 
Ap/p ~ 0.02 is the inlet pressure loss, the ratio of specific heats y= 1.4. Efficiencies are only 
valid for T T o > 80 kN 

HERRMANN, Steffen: Untersiichung des Einflusses der Motorenzahl auf die Wirtschaftlichkeit 
eines Verkehrsflugzeuges unter Beriicksichtigung eines optimalen Bypassverhaltnisses Berlin, 
Technical University, Institute for Aerospace Sciences, Department of Aircraft and Lightweight 
Design, Thesis, 2010 







































106 


Appendix B - Boeing Fuel Flow Method 2 


This method is adapted from Baughcum 1996 and Kim 2005. The aim is to generate a value 
for the emission index of NO x in function of the aircraft’s fuel flow. The method will be ex¬ 
plained in a step-wise manner. 

1. The aircraft’s engine should be looked up in the ICAO’s Aircraft Engines Emission Data¬ 
base (ICAO 2019). The ICAO values for fuel flow must be adjusted for installation effects on 
the aircraft. The correction factors are: 

• r — 1.010 for take-off 

• r — 1.013 for climbout 

• r — 1.020 for approach 

• r — 1.100 for idle 

2. In this step the adjusted fuel flow data and the corresponding EI_NO x values from ICAO 
2019 must be plotted in a log-log plot. An example is given: 



Figure B.1 : Example of the log-log plot of EI NOx in function of fuel flow (Kim 2005) 

3. The aircraft’s uncorrected fuel flow (in kg/s) must be calculated as follows: 

W f — f - V (5.1) 

With: V in km/s 

Then the corrected fuel flow must be calculated with the following equation from Boeing: 

Wff — (0 3,8 ■ e°- 2 MZ ) 


(5-2) 
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4. A linear trendline is plotted over the log-log plot in Excel. The equation from the trendline 
is then used to calculate the corresponding EI_NO x value and has the following general struc¬ 
ture: 


EInO x ,graph ' Wff T b 


(5-3) 


5. In this final step the calculated EI NOx graph value from step 4 is uncorrected in order to re¬ 
flect the at-altitude flight conditions. This is done by using different factors which take the ef¬ 
fect of humidity into account: 


k H 


/0.37318 ■ p v 

\p-0.6-p v 


0.0063 


(5.4) 


With: 


p v = 6895 ■ 0.014504 ■ 10^ 


(5.5) 


With: 


/? = 7.90298 ■ (1 - k T ) + 3.00571 + 5.02808 ■ log 10 (fc T ) 

+ 1.3817 ■ 10~ 7 ■ [l - I0 11344 (1 - fc r)] + 8.1328 ■ 10" 3 ■ [l0 3 ' 39149 ' (1 - fc7-) - l] (S. 6) 


With: 


373.16 
T + 0.01 


Finally, EI_NO x can be calculated: 


EInO x E lNO x ,g ra ph 



Q3.3 


(5.7) 
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Appendix C - Forcing Factor Data 


All values are compiled in Van Endert 2017, based on a graph from Schwarz 2011. In the 
original graph, data only started above 17 000 feet, but according to Schwarz the forcing fac¬ 
tor at lower altitudes can be set equal to the first available value above it. This principle is ap¬ 
plied to the data in red. 


Table C.1: Forcing factor for aircraft-induced clouds 


Altitude (ft) 

Forcing factor 

16000 

0,0284495 

17470,3 

0,0284495 

19547,9 

0 

21529,7 

0 

23511,4 

0,173542 

25525,1 

0,395448 

27506,8 

0,799431 

29456,6 

1,25178 

31598,2 

1,70982 

33547,9 

2,10526 

35529,7 

1,82077 

37543,4 

1,53343 

39557,1 

0,967283 

41538,8 

0,793741 


Table C.2: Forcing factor for short ozone 


Altitude (ft) 

Forcing factor 

16000 

0,469417 

17502,3 

0,469417 

19484 

0,55761 

21497,7 

0,620199 

23479,5 

0,711238 

25525,1 

0,711238 

27506,8 

0,813656 

29520,5 

0,930299 

31502,3 

1,00996 

33484 

1,13229 

35561,6 

1,42816 

37575,3 

1,62447 

39589,0 

1,8037 

41538,8 

1,93172 
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Table C.3: Forcing factor for long ozone and methane 


Altitude (ft) 

Forcing factor 

16000 

0,86771 

17470,3 

0,86771 

19484 

0,924609 

21497,7 

0,955903 

23543,4 

0,961593 

25525,1 

0,944523 

27538,8 

0,927454 

29520,5 

0,927454 

31534,2 

0,941679 

33516 

0,975818 

35561,6 

1,14083 

37543,4 

1,21479 

39589 

1,20341 

41570,8 

1,20341 
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Appendix D - Flight Time Estimation 


Flight time estimation is based on BAD A data from Eurocontrol 1998 and represented by the 
following four figures. 

TAS for climb and descent follows a pretty similar pattern for these aircraft (see figure D.l 
and figure D.2). The true airspeed will be approximated by two linear zones and a constant 
speed zone. A variation from 150 knots to 289 knots below FL100, a variation from about 360 
knots to cruise speed between FL120 and FL300 and a zone of constant cruise speed above 
FL300. 



A320 A340 B737 — B747 


Figure D.l: TAS during climb based on BADA data (Eurocontrol 1998) 



0 50 100 150 200 250 300 350 400 450 

Flight level 

— A320 —•— A340 B737 —*-B747 


Figure D.2: TAS during descent based on BADA data (Eurocontrol 1998) 


















Ill 


Rate of climb and rate of descent follow a more irregular pattern (see figure D.3 and figure 
D.4). A rough approximation will be made of an average ROC and ROD of 2 000 ft/min. 



A320 —•— A340 B737 —*-B747 - Approx. 


Figure D.3: Rate of climb based on BADA data (Eurocontrol 1998) 



—•— A320 — A340 —•— B737 —*—6747 - Approx. 


Figure D.4: Rate of descent based on BADA data (Eurocontrol 1998) 
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Based on this data, the following approximating equations are produced: 


Below FL100: 


TAS = 1.4 ■ FL + 150 


(D.l) 


At FL120: 


TAS = 360 


(D.2) 


Between FL120 and FL300: 


TAS = TAScr 360 ■ (FL - 120) + 360 
180 v ' 


(D.3) 


Above FL300: 


TAS = TAS cr 


(D.4) 


For the Airbus A320-200 this result in the average TAS during climb and descent as dis¬ 
played in Figure 6.6. 
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Appendix E - Case Study: Airbus A380-800 


In this appendix the procedure of calculating minimal fuel consumption, DOC, and environ¬ 
mental impact from Chapter 5 to 8 is repeated for the case of an Airbus A380-800. Only the 
final results for fuel consumption, DOC, environmental impact, and the combined model will 
be displayed. All aircraft data are found in Airbus 2005b. 


Mach number 



0,4 

0,45 

0,5 

0,55 

0,6 

0,65 

0,7 

0,75 

0,8 

0,85 

0,9 


3000 

8,0SE-05 

7,02E-05 

6.53E-05 

6,36E-05 

6.41E-05 

6,65E-05 

7.07E-05 

7.68E-05 

8.54E-05 

9,99E-05 

1.66E-04 


3500 

8,36E-05 

7.20E-05 

6,61E-05 

6,37E-05 

6.35E-05 

6,52E-05 

6,87E-05 

7.40E-05 

8.17E-05 

9.48E-05 

1.55E-04 


4000 

8.71E-05 

7.41E-05 

6,72E-05 

6,39E-05 

6.31E-05 

6.41E-05 

6.71E-05 

7.16E-05 

7.84E-05 

9.02E-05 

1.45E-04 


4500 

9,12E-05 

7.66E-05 

6,86E-05 

6,4SE-05 

6.30E-05 

6.34E-05 

6.57E-05 

6.95E-05 

7.55E-05 

8.61E-05 

1.37E-04 


5000 

9,57E-05 

7.95E-05 

7.03E-05 

6,54E-05 

6.31E-05 

6,29E-05 

6.46E-05 

6.78E-05 

7,30E-05 

8,25E-05 

1.29E-04 


5500 

1.01E-04 

8.28E-05 

7,24E-05 

6,66E-05 

6.36E-05 

6,28E-05 

6,38E-05 

6.64E-05 

7.09E-05 

7.93E-05 

1.22E-04 


6000 

1.06E-04 

8,66E-05 

7,49E-05 

6,81E-05 

6.44E-05 

6.29E-05 

6.33E-05 

6.53E-05 

6,91E-05 

7.66E-05 

1.16E-04 


6500 

1.13E-04 

9.08E-05 

7.78E-05 

7.00E-05 

6.55E-05 

6.34E-05 

6.32E-05 

6.45E-05 

6,77E-05 

7.42E-05 

1.10E-04 

JE 

7000 

1.20E-04 

9,57E-05 

8,12E-05 

7,23E-05 

6.70E-05 

6.42E-05 

6,33E-05 

6.41E-05 

6,66E-05 

7,23E-05 

1.05E-04 

O) 

7500 

1.28E-04 

l,01E-04 

8,50E-05 

7.51E-05 

6.89E-05 

6.53E-05 

6.39E-05 

6.40E-05 

6,59E-05 

7.08E-05 

1.01E-04 


8000 

1.36E-04 

1.07E-04 

8,94E-05 

7.82E-05 

7.11E-05 

6.68E-05 

6.47E-05 

6.43E-05 

6,55E-05 

6,97E-05 

9.70E-05 

< 

8500 

1.46E-04 

1.14E-04 

9.44E-05 

8,19E-05 

7.38E-05 

6,87E-05 

6.60E-05 

6.49E-05 

6.55E-05 

6,89E-05 

9.38E-05 


9000 

l,57E-04 

l,22E-04 

l,00E-04 

8,61E-05 

7.70E-05 

7.11E-05 

6.76E-05 

6.59E-05 

6,59E-05 

6,86E-05 

9.12E-05 


9500 

1.69E-04 

1.30E-04 

l,06E-04 

9,09E-05 

8.06E-05 

7,39E-05 

6.96E-05 

6.73E-05 

6,67E-05 

6,87E-05 

8.91E-05 


10000 

1.82E-04 

1.40E-04 

l,14E-04 

9.63E-05 

8.49E-05 

7.71E-05 

7,21E-05 

6.91E-05 

6,79E-05 

6,92E-05 

8.76E-05 


10500 

1.97E-04 

l,51E-04 

l,22E-04 

1.03E-04 

8.97E-05 

8.10E-05 

7.51E-05 

7.14E-05 

6,96E-05 

7,02E-05 

8.67E-05 


11000 

1.96E-04 

1.50E-04 

l,21E-04 

1.02E-04 

8.95E-05 

8,08E-05 

7.50E-05 

7.12E-05 

6,94E-05 

6,99E-05 

8.64E-05 


11500 

2,llE-04 

1.61E-04 

1.30E-04 

1.09E-04 

9.46E-05 

8,49E-05 

7,83E-05 

7.39E-05 

7.15E-05 

7.14E-05 

8.62E-05 


12000 

2,27E-04 

1.73E-04 

l,39E-04 

1.16E-04 

1.00E-04 

8,95E-05 

8,21E-05 

7.70E-05 

7.40E-05 

7,33E-05 

8.66E-05 


12500 

2.45E-04 

1.86E-04 

l,48E-04 

1.24E-04 

l,07E-04 

9.47E-05 

8,64E-05 

8.06E-05 

7.70E-05 

7.57E-05 

8.75E-05 

Figure E.1: 

Results fuel consumption A380 in kg fuel per NM, per kg of MTOW 




Mach number 



0,4 

0,45 

0,5 

0,55 

0,6 

0,65 

0,7 

0,75 

0,8 

0,85 

0,9 


3000 

0,174 

0,158 

0,146 

0,137 

0,130 

0,125 

0,121 

0,118 

0,117 

0,117 

0,116 


3500 

0,176 

0,160 

0,147 

0,138 

0,131 

0,125 

0,121 

0,118 

0,117 

0,117 

0,117 


4000 

0,178 

0,161 

0,148 

0,139 

0,131 

0,126 

0,121 

0,118 

0,116 

0,116 

0,117 


4500 

0,180 

0,162 

0,149 

0,140 

0,132 

0,126 

0,121 

0,118 

0,116 

0,115 

0,117 


5000 

0,182 

0,164 

0,150 

0,140 

0,132 

0,126 

0,121 

0,118 

0,115 

0,115 

0,117 


5500 

0,184 

0,165 

0,152 

0,141 

0,133 

0,127 

0,122 

0,118 

0,115 

0,114 

0,117 


6000 

0,185 

0,167 

0,153 

0,142 

0,134 

0,127 

0,122 

0,118 

0,115 

0,114 

0,117 


6500 

0,186 

0,169 

0,155 

0,143 

0,135 

0,128 

0,122 

0,118 

0,1 IS 

0,113 

0,117 

E 

7000 

0,188 

0,171 

0,156 

0,145 

0,136 

0,129 

0,123 

0,118 

0,115 

0,113 

0,117 

0) 

7500 

0,189 

0,173 

0,158 

0,146 

0,137 

0,129 

0,124 

0,119 

0,115 

0,113 

0,117 


8000 

0,190 

0,175 

0,160 

0,148 

0,138 

0,130 

0,124 

0,119 

0,116 

0,113 

0,117 

< 

8500 

0,192 

0,176 

0,162 

0,149 

0,139 

0,131 

0,125 

0,120 

0,116 

0,113 

0,117 


9000 

0,193 

0,177 

0,164 

0,151 

0,141 

0,133 

0,126 

0,121 

0,116 

0,114 

0,117 


9500 

0,195 

0,178 

0,165 

0,153 

0,143 

0,134 

0,128 

0,122 

0,118 

0,115 

0,117 


10000 

0,196 

0,180 

0,166 

0,155 

0,144 

0,136 

0,129 

0,123 

0,118 

0,115 

0,117 


10500 

0,198 

0,181 

0,168 

0,157 

0,146 

0,137 

0,130 

0,124 

0,119 

0,116 

0,117 


11000 

0,199 

0,182 

0,168 

0,158 

0,147 

0,138 

0,131 

0,125 

0,120 

0,116 

0,118 


11500 

0,199 

0,182 

0,169 

0,158 

0,148 

0,139 

0,132 

0,125 

0,120 

0,117 

0,117 


12000 

0,200 

0,183 

0,169 

0,159 

0,150 

0,140 

0,133 

0,126 

0,121 

0,117 

0,118 


12500 

0,201 

0,183 

0,170 

0,159 

0,150 

0,142 

0,134 

0,127 

0,122 

0,118 

0,118 


Figure E.2: Results seat-mile cost A380 in USD 
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Mach number 


0,4 

0,45 

0,5 

0,55 

0,6 

0,65 

0,7 

0,75 

0,8 

0,85 

0,9 

3000 

0,052 

0,022 

0,008 

0,003 

0,005 

0,011 

0,023 

0,039 

0,063 

0,102 

0,276 

3500 

0,060 

0,027 

0,011 

0,003 

0,003 

0,007 

0,017 

0,032 

0,053 

0,088 

0,248 

4000 

0,070 

0,033 

0,014 

0,004 

0,002 

0,004 

0,012 

0,025 

0,043 

0,075 

0,223 

4500 

0,082 

0,040 

0,018 

0,006 

0,001 

0,002 

0,009 

0,019 

0,035 

0,064 

0,200 

5000 

0,095 

0,049 

0,022 

0,008 

0,002 

0,001 

0,005 

0,014 

0,029 

0,054 

0,179 

5500 

0,110 

0,059 

0,029 

0,012 

0,004 

0,001 

0,004 

0,011 

0,024 

0,047 

0,164 

6000 

0,129 

0,072 

0,038 

0,018 

0,008 

0,003 

0,005 

0,010 

0,021 

0,043 

0,156 

6500 

0,150 

0,086 

0,048 

0,025 

0,012 

0,006 

0,006 

0,010 

0,019 

0,039 

0,146 

7000 

0,177 

0,105 

0,062 

0,036 

0,020 

0,012 

0,010 

0,012 

0,021 

0,039 

0,144 

7500 

0,204 

0,124 

0,075 

0,045 

0,027 

0,017 

0,013 

0,014 

0,020 

0,036 

0,134 

8000 

0,236 

0,147 

0,092 

0,058 

0,037 

0,024 

0,018 

0,018 

0,022 

0,037 

0,129 

8500 

0,280 

0,177 

0,115 

0,076 

0,051 

0,035 

0,028 

0,025 

0,028 

0,041 

0,132 

9000 

0,331 

0,213 

0,141 

0,096 

0,067 

0,049 

0,039 

0,034 

0,036 

0,047 

0,135 

9500 

0,385 

0,251 

0,169 

0,118 

0,084 

0,062 

0,050 

0,043 

0,043 

0,052 

0,135 

10000 

0,451 

0,297 

0,203 

0,144 

0,105 

0,080 

0,064 

0,055 

0,053 

0,061 

0,139 

10500 

0,540 

0,358 

0,248 

0,178 

0,133 

0,102 

0,083 

0,072 

0,067 

0,073 

0,149 

11000 

0,572 

0,379 

0,263 

0,190 

0,142 

0,111 

0,091 

0,079 

0,075 

0,082 

0,164 

11500 

0,685 

0,455 

0,319 

0,233 

0,177 

0,139 

0,115 

0,101 

0,094 

0,099 

0,181 

12000 

0,832 

0,554 

0,390 

0,287 

0,220 

0,176 

0,147 

0,129 

0,120 

0,123 

0,205 

12500 

1,000 

0,667 

0,471 

0,349 

0,270 

0,217 

0,182 

0,160 

0,148 

0,149 

0,230 


Figure E.3: Results environmental impact A380 per NM flown, per seat (normalized) 



Figure E.4: Results combined values A380 (normalized) 


The A380-800 shows a trend towards a higher Mach number and a slightly higher cruise alti¬ 
tude than the A320-200. The higher Mach number is because of the aircraft’s higher optimal 
cruise Mach number of Mach 0.85 compared to the A320’s Mach 0.78. The higher cruise alti¬ 
tude follows from the aircraft’s optimized design for higher altitude. The A380-800 is a wide- 
body long-range aircraft which is optimized to fly slightly faster and higher than a single aisle 
aircraft like the A320-200. The fact that it is built for long-range flight missions also explains 
why the seat-mile cost is significantly lower than the A320’s seat-mile cost. 



